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Chapter 2 deals with the theoretical model. According to 
this, a two-stage process sets in when an electric field is 
applied to the specimen. Firstly, Joule heating raises the 
sample temperature sufficiently to cause the metallic granules 
to melt. Subsequently, the electric field stretches the molten 
granules against the surface tension force and at the threshold 
field, the granules touch each other thereby causing the material 
to switch to the "on" state. A one dimensional heat equation has 
been solved numerically to find the electrical field necessary to 
raise the temperature to the effective melting point of bismuth 
particles. However, the critical field necessary to stretch the 
molten granules so that they touch each other is found to be 
larger than the previous one and hence determines the overall 
switching field. The values of the latter for different particle: 
diameters and separation to diameter ratios have been calculated-! 
These compare reasonably well with experimental data obtained in ^ 
the cases of silicate and vanadium phosphate glasses containing 
bismuth particles. 

In Chapter 3, the voltage -current characteristics of 
vanadium phosphate glasses containing bismuth granules especially 
in the negative resistance region have been computed using the 
model developed in Chapter 2. It is shown that according to the 
mechanism proposed first a conducting region is formed at the 
specimen centre due to the metal spheres touching each other. A 
regenerative process then sets in as a result of increase in the 
current flow through the specimen ultimately leading to the "on" 
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state when the percolation of the metal phase in the glass 
matrix is attained. One dimensional time dependent heat 
equation has been solved numerically at different instants of 
time to delineate the corresponding temperature profile within 
the sample and the concommitant changes in its resistance and 
the current flow through it. Experimental measurements on well- 
characterized speciuBns in the have been 

carried out and the results are compared with the theoretical 
values. Satisfactory agreement between the two has been obtained 
It is observed that increase in sample thickness leads to a 
threshold switching. Eor a memory switched sample measurement 
of its resistance in the "on" state as a function of temperature i 
shows that in the range 75°C - 110°C it reverts back to the off 
state through a voltage controlled negative resistance in the 
range 0.2 to 0.5 volt. All these have been explained on the 
basis of the present theory. I 

In Chapter 4, the kinetics of switching has been investi- 
gated by solving numerically fluid flow equations corresponding j 
to the movement of boundaries of the molten granules during the j 
stretching process. Switching time has been measured experi- 
mentally by applying electrical pulses on samples in ^ 

Bi 20 ^ system at room temperature and monitoring changes in their j 
resistance. Experimental values (in the range 50 - 80 micro- 
seconds) are found to be orders of magnitude higher than the i 

theoretical prediction ( ^ 10”^ second). However, the ; 
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switching time computed from the solution of time dependent 
heat equation is in reasonable agreement with experimental data. 

It is concluded, therefore, that the kinetics is controlled 
primarily by the Joule heating step. 

In Chapter 5 » the optical properties of vanadium phosphate 
as well as sodium borosilicate glasses containing bismuth 
granules of dimensions of the order of a few hundred angstroms 
have been investigated. The micro structures of samples having 
different compositions and subjected to various reduction treat- 
ments have been characterized by transmission electron microscopy. 
Optical absorption in the visible range has been measured using 
Carry 1 7D Spectrophotometer. Vanadium phosphate glasses con- 
taining bismuth show an absorption peak around (450 +10 nm) 
while the silicate glasses containing bismuth show two peaks in 
the ranges 500 nm - 530 nm and 420 nm - 430 nm respectively. The I 
peak position in V 20 ^-P 20 ^-Bi 20 ^ glasses is predicted in fair j 
agreement with experiment by Maxwell-Garnet t (MG), texwell- 
Garnett as extended by Polder and Van Santen (MG-PVS) and i 

j. 

Bruggeman (BE) effective medium theories with MG theory givii^ I 
the best fit. In the case of silicate glasses containing bis- | 

I 

muth the predicted positions do not agree so well with experi- 
mental results. The possibility of a fraction of the metal 
spheres being present in the form of interconnected metallic 
chains and/or f.c.c. clusters in the present microcomposites 
has also been explored. A typical computation has been carried 
out by using a constrained optimisation technique with some I 



assumed fractions of spheres, f.c.c. clusters and single 
strand chains as the optimising parameters. It is observed that 
in a vanadium phosphate glass containing bismuth quantitative 
agreement between the computed and experimental a values within 
an error of 4% over the entire spectral range can be achieved 
if the metallic phase is assumed to consist of 26% spheres 
0.6% single strand chain and 73 *6% f.c.c. clusters. 

In Chapter 6, the main conclusions drawn on the basis of 
the entire investigation are summarised. 



CHAPTER 1 


IHTRODUGTIOH 

1 « 1 Applications of Class-Metal Microcomposites 

Glass-metal microcomposites consist of a glass matrix 
with a distribution of metallic particles in the range of a few 
hundred angstroms in size. The metallic particles may b© preci- 
pitated either during melting, controlled cooling of the melt 
or during subsequent heattreatment depending on the glass com- 
positions and melting conditions. The applications of glass- 
metal microcomposites may be broadly classified as follows^ . 

1 . Photosensitive glasses 

2* Photosensitively nucleated glass-ceramics 

3. Photochromic glasses 

4. Polychromatic glasses 

5. Photothermal conversion (cermet) 

6. Solar control coatings (electro float process) 

Some of the other possible applications will be : 

1 . Electro-conducting glasses and fibres 

2. Memory switching 

3. High strength glasses and fibres 

A brief review of the present state of research of the glass- 
metal microcomposites is presented in the following sections. 
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1 .2 Photosensi t ive G- lasses 

In photosensitive glasses, metallic granules of colloidal 
dimensions are precipitated by subjecting the glass to certain 
heattreatment schedules. The precipitation is greatly enhanced 

if the glasses are exposed to UV radiation prior to heat- 

2 

treatment . 

3 

Weyl considers that copper, silver and gold may be 
present in glass in tme solution and they get precipitated 
either by controlled cooling of the molten glass or by subse- 
quent heating of the same. The relative proportions of ions 
(Gu^, Ag"^ and Au^"^) and the corresponding metallic atoms 
depend on the melting condition and the composition of the 
glass. When thb molten glass is quickly cooled the high 
temperature state is retained and the glass remains colourless. ; 
The glass- strikes the colour when it is subjected to a heat- 
treatment. The mechanisms involved are as follows. During 
heattreatment the reducing agents, e.g., SnO, As20^,Sb20^ etc- 
get oxidised after reducing suitable ions in the melt to their 
metallic counterparts. These metal atoms along \jith those 
which are already present in the matrix nucleate and grow to 
colloidal dimensions. These metallic granules are responsible 
for the selective absorption of light and so colour is developed! 
In this connection, it should be mentioned that platinum in the ! 
form of colloidal dispersion has not been used as colourant, ! 
because this leads to greyish tint. However, the importance of 
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colloidal platinum is that it can he used as a nucleation 
catalyst to prepare glass-ceramics. In contrast to copper, 
silver and gold no heattreatment is necessary to form colloidal 
platinum^. The platinic chloride which is usually added as 
the source of platinum breaks during melting to form platinum 
colloids in the melt. 

The earliest photosensitive glasses contain copper as the 
photosensitive element. In later development it has been shown 
that silver and gold can also act as photosensitive elements. 

A range of photosensitive glass compositions containing these 

5 

elements have been worked out by Stookey . The fundamental 
mechanism of the photosensitive process involves two steps. 

The glass is first irradiated selectively by placing a photo- 
graphic negative to form the latent image. This is achieved by 
reducing the ions to their atomic states as follows : 

Cu"^ + e-»'Cu° (l.l) 

The electron necessary for this reduction reaction is obtained 

-L p4- 

by photo-assisted oxidation of Cu to Cu ions, according to 
the following reaction, 

Cu"^ + h V Cu^'*' + e (1»2) 

The irradiated glass is then subjected to a heattreatment 
in the temperature range of 500°G-600°C for a suitable period 
of time. Since the irradiated portion of glass contains more 
metallic atoms the metallic phase nucleates and grows to give 
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colour and thus the latent image is developed. 

It has been shown by Stookey^ that inclusion of cerium 
oxide in the glass batch greatly improves the photosensitivity. 
This is essentially due to the ease with which electrons are 
made available by the reaction, 

Ce^^ + h Ce^'*' + e (1 .3) 

G-old and silver are introduced in the photosensitive glass 
batches as their chlorides. The glass batches are melted in 
oxidising condition. The amount of silver present is equivalent 
to 0,05% of 0.3% AgCl, while the concentration of gold varies 
in the range of 0.01% to 0.1%. On the other hand glasses con- 
taining 0.05 to 1 ,0fo of copper (calculated as GU 2 O) as photo- 
sensitive element are melted in reducing atmosphere. The 
glasses contain small amount of Sb 20 -^ (about 0.1%) along with 

cerium oxide. The beneficial effects of Sb20^ lie in the fact 

5+ 

that it provides electrons by being oxidised to Sb ionic 
state . i 

*3 P hotosensitively Nucleated Glass-ceramiCs 

The process of photosensitively nucleated glass -ceramics 
involves three steps. Firstly, the chilled glass is exposed to 
U¥ radiation so that the fomation of metal atoms takes place. 
Secondly, the glass is heattreated to a temperature which depends 
on the system used, for a certain period of time so that sub- 
microscopic metallic particles get precipitated by a process of 
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nucleation and growth. Lastly, the glass is further heated 
to higher temperatures when the metallic granules act as 
heterogeneous sites for controlled crystallisation of glass. 

This is one of the early methods of producing glass-ceramics. 

7 

Stookey has studied many glasses of the Li20-Al20^-Si02 
system. These glasses contain either gold in the range of 
O.OOlfo to 0.003f° or copper in the range of 0.001% to 1.0% 
(computed as CU 2 O) or 0.001% to 0.5% of silver (calculated as 
AgCl) . The glasses are exposed to UV radiation followed by 
heattreatment at 500°-540°C and subsequently heattreated at 
temperatures around 950°C to produce the glass-ceramics. These 
materials have good mechanical strength and electrical proper- 
ties. In a later development it has been shown by Stookey 

8 o 

et al. that there is a minimum size of metal namely 80 A below ^ 

which they cannot act as nucleation catalysts to crystallise 

4 

lithium metasilicate. Rind one has studied the use of metallici 

platinum as heterogeneous nucleation catalyst to crystallise 

lithium silicate glasses. It has been demonstrated that an 

optimum amount of, platinum (0.004% to 0.007% at 600°G) is 

needed to get maximum crystallisation. It has been shown that 

there is an optimum radiation dose necessary for the photosensi-i 

tively nucleated gold containing lithium aluminosilicate glass-; 

2 

ceramics to attain the maximum mechanical strength . ; 

1.4 Photoehromic Glasses 


¥hen photochromic glasses are exposed to sunlight they 
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darken and fade when the light is removed. This is a rever- 
sible process. A similar process called solarisation has been 
known for many years. In contrast to photochromic glasses the 
latter are sensitive to high energy electromagnetic radiation. 

9 

Recently, Araujo has reviewed different aspects of photo- 
chromic glasses. 

1 0 

Armistead et al, have developed alkali alumino- 
borosilicate system which has perhaps been . . most thoroughly 

studied. The other systems investigated are alkali and alkali 

11 12 
earth borates , heavy metal aluminoborate s to name a few. 

The silver halide crystal which gets precipitated, induces 

photochromic properties. The mechanism involved is believed 

to be quite analogous to that operative in photographic films. 

When the photochromic glass is irradiated with light silver 

halide is decomposed into silver and halogen. Unlike the 

photographic film, the liberated halogen is trapped into close 

proximity of silver due to the rigidity and impermeability of 

glass. When light is removed the silver and the halogen atoms 

recombine and revert to transparent colourless silver halide 

crystals. The reversible reaction, which causes darkening and 

1 5 

fading of photochromic glasses can be represented as follows 

'^''1 

AgX * Ag° + X° (1 .4) 

hv^ 

In many photochromic glasses, the absorption band is found near 

14 

500 nm and this is attributed to colloidal silver by Moriya . 
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Seward suggests that ellipsoidal silver particles are res- 
ponsible for optical absorption in the visible range-. The 
theoretical model explaining the behaviour of photochromic 

Q 

glasses can be found in a review article by Araujo . 

If the glass is cooled very rapidly from molten state, 
photochromic property is not developed because no silver halide 
Particles are precipitated. To develop photochromism a sub- 
sequent heattreatment is necessary. There is an optimum range 
of silver halide particle size, namely 80-150 2 ., If the 
particles are smaller than the lower limit the glasses do not 
darken at all. On the other hand if the particles are larger 
than the upper limit the glasses become hazy. Since the 
increase of subsequent heattreatment temperature increases the 
particle size usually silver halide containing photochromic 
glasses are heattreated in the temperature range of 400°-800°C 

9 

for 1 5 minutes to 4 hours . After the precipitation of silver 

halide particles, further low temperature heattreatments have 

1 6 1 T 

profound influence on photochromic properties ’ . The combined 

influence of compositions and different heattreatment tempera- 
ture for different periods of time of some glasses in the alkali 

1 8 

aluminoborosilicate system has been studied by Hares et al. . 

It indicates the importance of base glass compositions and 

thermal history on the darkening and fading characteristics of 

19 ' 

the photochromic glasses. Gliemeroth has discussed the 
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properties of photochromic glasses as a function of temperature 

and radiation dose. This is important from the point of view 

of applications in different weather conditions. The presence 

of copper in small amount ( 0.02%) increases the darkening 

1 3 

sensitivity very effectively as it augments the^ production 
of neutral silver atoms according to the reaction, 

hv^ 

Ag^ + Cu"*" t- Ag° + Cu^'*’ (1.5) 


hv 


2 


20 

Similarly the presence of cadmium only in presence of copper 
enhances the darkening sensitivity. Along with the composition 
as mentioned earlier, the melting conditions can affect the 
photochromic properties as they strongly influence the retention 
of halide and oxidation states of copper and silver in the glass 
me It . 


We have, so far, discussed the photochromic glasses con- 
taining the silver halide particles. However, highly reduced 

21 22 

homogeneous alkali silicate glasses , Gd 0 -B 20 ^“Si 02 glasses , 
Ha 2 O -Al^O^ .23102 (S being chloride or sulfide) glasses 

show photochromic behaviour also. Glasses containing suspension 

„ 24 

of crystallites-hackmanite , silver molybdate , thallium 

20 

chloride ^ and copper cadmium halides show photo chromism. 

The thermally darkenable photochromic glass (TDPC) is 
quite different from the usual photochromic glasses. Unlike 
the latter the TDPC glasses are dark at room temperature and 
become darker with further heating. When exposed to visible 
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light the glasses bleach to colourless state at a rate 
dependent on the intensity of light. On removing the source of 
light the TDPC glasses returns to darkened state at a temperature 
dependent rate. Seward III has described the TDPC proper- 
ties of lanthanum borate based glasses containing copper doped 
silver halide crystals. The composition ranges of these 

t 

lanthanum borate glasses can be found in the patent 
27 

literature 

26 

Seward III has observed high temperature bistable optical 
states of TDPC glasses. When the glasses are heated in the 
temperature range of 400°-450°C, they are bleached to a colour- 
less state. When the temperature is lowered to the range of 
300°-350°C, they spontaneously go to a darkened state and 
do. not bleach despite the presence of intense light. These 
two bistable states are believed to arise due to melting and 
precipitation of silver halide crystals^S. 

The proposed mechanism of TDPC glasses is the precipita- 
tion of colloidal silver from the molten silver halide droplet 
of nonstoichiometf ic composition during the crystallisation of 
the latter. Copper may act as either Gu"^ hole trap or Cu2+ 
electron trap and may take part in the optical bleaching or the 
thermal darkening. On the other hand its only role may be the 
controlling of oxidation states and hence producing more neutral 
silver atoms required for the precipitation of colloidal silver 
from molten halide droplets. 
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^ • 5 Polychromatic Glasses 

The polychromatic glasses are a special kind of photo- 
sensitive glasses. A variety of colours can be generated by 
the application of controlled thermo-optical treatments. 

The sequence of steps involved to form the polychromatic 
29 

glass is as follows . Firstly, the glass is exposed to UV 

radiation to precipitate the small silver particles and this 

3 + 

is assisted by the presence of Ce photosensitisers . In the 

second step the glass is subjected to a heattreatment such that 

the silver particles act as heterogeneous sites for the nuclea- 

tion and growth of sodium halide microcrystals. It has been 

30 

shown by Stookey that the sodiim halide crystals have cubic 
base with different degrees of pyramidal protrusion. A second 
UV exposure followed by heattreatment serves to produce silver 
at the tips of the sodium halide crystals. The amount of 

'T /-V *7 A 

silver is related to the type of colour generation^ * . 

The colour which is obtained is primarily a function of the 
first UV exposure . A strong UV exposure produces a large number 
of silver particles and hence a large number of sodium halide 
crystals is precipitated as a result. The second UV exposure 
produces a relatively fixed amount of silver. This amount of 
silver is then distributed over the tips of the crystals. So, 
if the number of sodium halide crystals is less it produces 
longer and larger specks (higher aspect ratio) compared to the 
case when the number of crystals is' large. This aspect ratio 



controls the generation of colour in polychromatic glasses. 

1 .6 Photothermal Conversion 

The efficiency of photothermal conversion of selective 
black absorbers can be increased by increasing the absorptivity 
over the solar radiation spectrum and by decreasing the 
infrared emissivity. In principle this can be obtained in a 
single material. However most of the selective black absorbers 
now available is produced by depositing a thin film of metal- 
dielectric composite which is transparent in the infrared range 
but highly absorbing in the visible range, on a metal substrate 
which has low infrared emissivity. The commonly used photo- 
thermal conversion materials are electroplated Ni-black, Cu-black 
on metal sheets. The operating temperature required for photo- 
thermal conversion may vary from 400°G - 1000°G. However, the 
Hi-black and Gu-black are not stable above 200°G while Gr-black 
is not much stable above 200°G. The MgO/Au^ cermet film has 
been produced by rf sputtering on metallic sheet. It has 
excellent solar absorptivity viz., over 0.9 and its infrared 
emissivity is less than 0.1 . With Mo metal sheet the selective 
black absorber is stable upto 400°G. The Gr20^/Gr^^ rf 
sputtered film on metal has solar absorptivity and infrared 
emissivity similar to I’'igO/Au cermet film. Its behaviour remains 
essentially unchanged even after heating at 300°C for 60 hours. 
The co-sputtered Hi/Al20^ cermet film on metal has also 
similar absorptivity and emissivity but it is thermally more 
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stable (stable upto 500°C) as compared to other films. 

3 S 

Recently, Riklasson et al. have shown by detailed optical 
study on evaporated nickel particles that these are stable 
upto 600°C in air without lossing their selectivity. This 
makes the use of nickel grains very useful for photothermal 
conversion of solar energy. 

1.7 Electrofloat Glasses 

Electrofloat process involves the electric field induced 
ion-exchange of alkali containing glasses, one of the ion- 
exchange species being copper. The ion-exchanged glass is then 
subjected to a reduction treatment to precipitate the metallic 
particles. This process has been developed by Pilkingtonsto 
produce the solar control coating. 

^ • S E lectro-conducting Glasses 

1 ,8..1 Metallic regime 

A number of methods has been developed to obtain a highly 

conducting layer on or within the surface of the glass. These 

can be divided under two broad catagories. In the first, the 

layer is simply applied on to the surface of the glass. 

Platinum, gold^ or copper salts are heated to reduction on 

37 

the surface, metals are vaporised on to the surface , silver 
oxides and metal powders are fired on with the help of 
auxiliary fluxes and lacquers containing suspensions of flaky 
metals^ are applied. The second category involves techniques 
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whereby various glasses containing lead oxide and other metal 
oxides are reduced in hydrogen gas to induce a conducting 
layer within the glass surface^^ 

Some alkali containing silicate glasses develop hi^ 

41 42 

surface conductivity ’ when they are subjected to an ion- 
exchange treatment followed by reduction. The process involves 
the dipping of the glass in molten silver nitrate at about 
300°C for certain period of time so that alkali silver ion- 
exchange takes place. Then the glass is reduced in hydrogen 
gas in the temperature range of 250°C to 400 °G so that silver 
metal atoms nucleate and grow to form the interconnected chains 
of metal and hence a highly conducting layer is induced within 
the surface of the glass. This is confirmed by electron micro- 
structural analysis. The surface resistance is typically in 
the range of 0.15 u /p to 4*1 ^ /o depending on temperature 
and diiration of the reduction treatment. It is found that 
prior to the thermo-chemical treatment it is necessary to make 
the glass surface sufficiently rough by grinding with 1 20 mesh 
silicon carbide powder to induce highly conducting layer. It is 
believed that the rough surface provides the potential sites for 

silver metal to nucleate and grow during the reduction treatment. 

43 

It has been shown recently that the ceramised GaO-S 102-159-2*^" 
Bi 20 ^ systems develop high surface conductivity when they are 
ion-exchanged in molten silver nitrate followed by reduction in 
hydrogen gas. The typical surface resistance varies in the range 



of 0.08 y/a to 14.76 Q /o depending on the compositions and 
reduction temperature and time. The glass-crystal interface 
of the partially crystallised glass provides the potential 
sites for the silver metals to nucleate and grow during 
reduction treatment. Hence, surface roughness is not necessary 
in this process of generating highly conducting layer. 

1 . 8.2 Dielectric regime 

When oxide glass systems containing metals are melted 
under ordinary or reducing conditions, they develop a micro- 
structure very similar to that of cermet films produced by 
sputtering or evaporation. Mixtures of oxides in the 2^3 

system containing Ag2^ or , Na20-B20^-Si02 system 

containing Bi20^ or pure selenium"^^ and Pb0-B20^ system^^ 
containing Bi20^ when melted and cooled rapidly induce preci- 
pitation of metallic granules of silver, bismuth or selenium 
in the corresponding glass matrix. The precipitated phases 
have been confirmed by selected area electron diffraction. The 
particle size found by electron microstructural analysis is of 
the order of a few hundred angstroms. Compared to base glasses 
the metal containing glasses show a decrease in resistivity. 
These glasses containing a distribution of metallic particles 
show semiconducting behaviour. Prom the temperature dependent 
conductivity measurements it has been concluded that the 

mechanism involved is the tunnelling of electrons between 

45 47 

conducting metallic islands . 
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The borosilicate system^ containing Bi20^ and/or Ag^O 

have been melted in a reducing condition. Fibres have been 

49 

drawn from single hole bushing made of pure Al 20 ^ , at a fiberi- 

sing temperature of around 1100°C, The fibres are found to have 
a microstructure consisting of bismuth and/or silver metallic 
particles of sizes ranging from 100 to 500 dispersed in the 
glass phase. The conductivity is electronic in nature and the 
mechanism is found to be electron tunnelling between the condu- 
cting metallic species similar to what has been observed in the 
case of bulk glass-metal microcomposites. 

1 . 9 Memory Switching,' 

In certain amorphous systems, when an applied electric field 

to the sample crosses a certain critical value the sample goes 

from a high resistance state to a low resistance state . This 

phenomenon is known as switching. There are two types of 

switching - threshold and memory. When the applied electric 

field is removed the threshold switch reverts back to the off- 

state. On the other hand the memory switch remains in the on- 

state even after the removal of the electric field. Alkali 

borosilicate glasses containing bismuth oxide when subjected to 

ion-exchange in molten silver nitrate followed by a reduction in 

hydrogen are found to have micro structures consisting of meta- 

50 

llic particles distributed in a glassy matrix. This system 
is found to exhibit memory switching with a reduction of resi- 
stance of six orders of magnitude when the applied electric field 
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crosses the critical value of 30 V/cm. When a 5-”'volt pulse of 

10 msec, duration is applied to the sample it goes from the 

off to the on-state. However, on application of a 50-volt pulse 

of 1 0 ju, sec. duration the sample switches back to the off-state. 

-5 

These glasses have been evaporated in a vacuum of 10 torr on 

to alumina substrates to form thin films of a few micron thick- 

51 

ness. Such films also show memory switching. When the sample 
in the on-.state is heated above a temperature (around the melting 
point of bismuth metal) it reverts from the on to the off -state , 
Recently, V^O^-P^O^ glasses^^ containing metallic bismuth or 
silver have been reported to show memory switching. The memory 
Switching action has been explained on the basis of a particle 
stretching model, 

1.10 High Strength Glasses and Fibres 

48 

Alkali borosilicate glass fibres containing bisnuth 

and/or silver are found to have higher strength compared to those 

drawn from the base glass composition. However, the improvement 

is small. This increase in strength is believed to arise due to 

the deflection of cracks along the glass-metal interface during 

52 

the course of crack propagation , All fibres containing metallic 

silver particles are found to have higher Young's moduli! compared 

to the base glass fibres. This arises due to a higher modulus 

53 

of silver. Si02“B2^3"'^2®~^^2*^3*"'^®2^3 glasses containing silver 
are found to have higher rupture moduli! as obtained from the 
three point bending measurement compared to base glasses. 
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Similarly alkali aluminoborosilicate glasses containing 

\ 

bismuth have higher rupture modulii and Young's modulii compared 
to base glasses. 

^ ^ 1 Present Work 

Memory switching observed in silicate and vanadium phos- 
phate glasses containing metallic granules of bismuth has been 
qualitatively explained as due to the stretching of these 
particles under the action of an electric field. The objectives 
of the present thesis are the following : 

1 . To develop a detailed model to explain quantitatively 

the memory action in oxide glasses containing microgranules 
of metallic bismuth, specifically to calculate the threshold 
field as a function of different microstructural parameters. 

2. To compute the current-voltage characteristics of these 
systems using this model especially in the negative resi- 
stance region, 

3. To estimate the switching time for different micro- 
structural configurations. 

4. To carry out suitable experimental measurements to test the 
validity of the model. 

5(a) To undertake detailed optical characterization of these 
glass-metal composites. ■ 

(b) To apply the different models of inhomogeneous media to 
explain the optical absorption data. 



CHAPTER 2 


* 

PARTICLE STRETCHING MODEL AND THRESHOLD FIELD IN 
YMADim PHOSPHATE AND SILICATE GLASSES 
CONTAINING HISMUTH GRANULES 

2 . 1 Introduction 

Preliminary reports on switching in chalcogenide glasses 

ET cr jr /ET 

first appeared in the literature nearly two decades ago^^’^ . 

57 

Ovshinsky’s paper in 1968 described the prototype of a relia™ 

ble potential device and subsequently extensive theoretical 

and experimental investigation have been undertaken in various 

laboratories to understand the phenomenon. The switch described 

by Ovshinsky is a sputtered film of a typical thickness of 

about 1 urn sandwitched between two metal or graphite electrodes. 

The material used is an amorphous semiconductor with an activa- 

58 

tion energy '^0.5 eV and optical’ band gap '^1.0 eV . The 
physical properties of these materials have since been widely 
discussed in literature^^"'^^ . 

Switching can be of two types •" threshold and memory. 
Figures 2.1(a) and 2.1(b) show the schematic representations of 
I-"V characteristics of a threshold and a memory switch res- 
pectively. Both of these switches go from a high resistance 
state (off.^state) to a low resistance state (on-state) via 

A paper based on this work has been published in J.Appl.Phys. 
51 , 3896 (1980). 





(a ) 



Fig. 2.1 Schematjc representations of I-V charaGt- 
eristics of (a) threshold (b) memory switch 
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negative resistance when the applied electric field crosses a 
critical value. This is denoted as in figure 2.1. In 
contrast to the threshold switch, a memory switch remains in 
the on-state even after the removal of the applied field. 

However, a minimum holding current is required to retain 
the on-state of the threshold switch. When the applied field 
is removed, the threshold switch returns to the off-state. A 
strong electrical pulse with short duration is required to “bring 
the memory switch from the on-state to the off -state. 

Whether a chalcogenide glass exhibits threshold or memory 

67 

switching is primarily determined by its composition . .A 

typical- threshold composition is As^gTe^gSi^ 2^® ^0 so-called 

57 

STAG- glass) . Usually most of the memory compositions are 
close to Ge^^TCg^X^ (X being As,S,Sb etc. or a combination^^ ^ , 
The composition dependent switching characteristics of chal- 
cogenide glasses can be understood qualitatively from their 
glass forming ability. The glass forming regions in the two 
ternary systems Ge-Te-As and Si-Te~As have been experimentally 
determined'^'^ ^ . Prom these Owen et al.^*^ have delineated the 
glass forming region of four component system, Ge-As-Te-Si by 
interpolation. It is reported that the so-called STAG glass 
falls at the centre of the glass forming region and hence the 
STAG glasses are resistant to crystallisation. On the otherhand 
the memory glass composition Ge^^Tsg^As^ falls near the edge 
of the glass forming boundary. So they easily get crystallised 
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above the glass transformation temperature (fg) ^^.t at the 
same time can be quenched quickly from molten state to form 
glass. Thus the choice of memory composition should be a 
balance between the ease of crystallisation above T^ and 
tendency to glass formation. The threshold composition 
should be such that it makes the glassy state as stable as 
possible. 

Two types of theories - thermal and electronic have 
been proposed to explain the switching behaviour. The classi- 
fication depends on whether the on-state is achieved by heating 
(thermal model) or by non-equilibrium enhancement of carriers 
(electronic model). In the thermal model, the application of 
voltage across the sample causes the temperature to rise to a 
critical value due to Joule heating. After a certain time when 
the critical temperature is exceeded no stable temperature .. 
distribution is possible and any fluctuation causes the tem- 
perature to rise very fast until an upper limit, set by the 

72 

external circuit . resistance, is reached . The result is the 

formation of the hot conducting channel in the on-state and the 

associated negative resistance state prior to switching. 

Earlier investigations on thermal instability have been reported 

73 74 

by Lueder et al, and Ridley . Recently, many 

rrj no 7 -wRO 

authors” ' ’ ' have applied the thermal model to explain 

the threshold switching characteristics. Some author 
have also incorporated the field dependent conductivity along 
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•with the temperature dependence to explain the switching 
characteristics in chalcogenide glasses. 

According to the thermal model, the system returns to the 
off-state when the slope of I-Y characteristics of the sample 
is equal to that determined by the load line of the circuit. 
Therefore, depending on the value of the external resistance 
the cut-off points should vary. Indeed this is observed by 

O 4 

Hughes et al. . The sample studied is a thick film. On the 

QO 

other hand Esquda et al. have observed no matching of slopes 
for thin samples. A system where switching is possibly caused 
by the thermal mechanism is the vanadate glass system studied 
by Higgins et al, , The switching in these glasses is attri- 
buted to a metal-insulator transition at 68°C. 

84 

The electronic model was first proposed by Mott and 
85 

Henisch . The theory states that in the on-state s'ufficient 

double injection of carriers takes pla.ce to ensure the presence 

of degenerate electron and hole gas. The carrier density is 

10 /c.c. and the Permi energy of both carriers lies on the 

86 

extended side of the respective mobility edge . Thus the 
potential drop across the sample is small and obviously the 
holding voltage is slightly greater than the mobility gap. A 
number of authors^"^””^^^ have used the electronic model to 
explain the chalcogenide switching characteristics. 

To maintain the high density of electron and hole gas in 
the on-state, the supply of carriers should be there to 
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balance the loss of carriers by recombination and by diffusion 
out of the glass sample. Thus the minimum holding current 
represents the condition when the replacement of carriers fails 
to balance the loss and so the switch changes from the on-state 
to the off-state^®. 

In memory switching the on-state is induced by the precipi- 
tation of the conducting crystals in the glass system. Ovshinsky 
93 

et al. have shown that strong illumination can induce a similar 

94 

behaviour. Pritzsche et al. have observed by calorimetric 

measurement that there exists an exothermic peak between 

20d^-300'^C and an endothermic peak at higher temperature for 

chalcogenide glasses showing memory switching. Therefore it is 

hot entirely clear whether the crystallisation is induced by 

86 

heating or by a high density of carrier injections . To bring 
the memory switch from the on-state to the off -state a high 
amplitude pulse with short duration is applied such that the 
crystals get melted and quenched to the glassy state. 

As described in Section 1.9 of Chapter 1, alkali borosilicate 

glasses containing bismuth show memory switching after ion- 
50 

exchange and reduction. It has recently been reported that 

44 

glasses containing bismuth and silver and the 

evaporated thin films of alkali borosilicate glass containing 
51 

bismuth show memory switching. The typical switching field 

4 / 

is of the order of 10 7/ cm. Unlike the chalcogenide glasses 
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where memory switching is explained by ciystallisation of 

glassy phase, the memory switching in glass-^netal composite 
44 

systems has been explained qualitatively on the basis of 
particle stretching model. The model essentially states that 
if the metallic particles are in molten condition, the applica- 
tion of an electric field causes the particles to stretch which 
is being .opposed by a surface tension force. Thus a critical 
electrical field causes the particles to stretch and form inter- 
connected chains thereby forming a low-resistance on-state. So, 
in this chapter we shall develop this model to predict quanti- 
tatively the critical field for switching in the glass Hnetal 
micro -particulate system mentioned above. 

2.2 Theoretical Model 

When an electric field E is applied to a dielectric matrix 

containing uniformly distributed metallic particles, the field 

/ 

induces opposite charges on the surface of the conducting 
species. If these metallic particles are in a molten state, 
two forces acting opposite to each other are operative. One is 
the coulombic attraction arising between two oppositely charged 
neighbouring metallic grains. This attraction force tends to 
stretch the molten metallic granules. The other is the surface 
tension force which opposes the stretching force . Therefore, 
depending on the magnitude of the applied field an equilibrium 
will be reached when the two forces balance each other. The 
threshold field then corresponds to the situation when two 
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neighbouring grains touch each other. ¥e assume that the 
metallic particles are spherical in shape in accordance with 
the experimental observation. Further we make the following 
assumptions for the simplicity of calculation : 

1 ) The electric field is uniform throughout the matrix and 
is normal to the electrodes. 

2) When the molten sphere is stretched it becomes ellipsoidal 
with major axis forming the axis of revolution. 

2.2,1 Determination of electrostatic stretching force. F ^^ 

To find the expression for the electrostatic force (^’^2.^ 
between the two polarised ellipsoids due to the presence of an 
electric field E, the left hand ellipsoid is displaced by 6x 
in the negative direction of x-axis as shown in figure 2.2. 

The work done against F^^ is given by 

6w = Fg^ 6x . (2.1) 

This work is equal to the electrostatic energy stored in the 
volume (6v) swept by the ellipsoid. Since electrostatic energy 

p 

density is given by , 

F^^6x = -jee^E^dv (2.2) 

where e and are the relative dielectric permittivity of the 
glass-metal composite and the dielectric permittivity of free- 
space respectively. 




27 


2. 2. 1.1 Determination of volume 5v 

let us consider an elemental circular strip of area dA 
at the point A in figure 2.2, The area dA is given by 


dA 



1 + y'2 dx . 


(2.3) 


Due to the movement of strip by 6x the volume (dv) is given by 


dv = dA n*i 6x 


(2.4) 


where n and i are the unit normal vector and imit vector along 
the X axis respectively. 


Dor an ellipsoid with major axis as the axis of revolution S 
is given by 


n 


= (I 


X 


a 


•f 




(2.5) 


where, a and b are respectively the semi -major and semi-minor 
axes of the ellipse respectively. 


Substituting dA and n respectively from equations (2,3) and 
(2.5) in equation (2.4) and then integrating the resulting 
expression from x = o to x = a, we get 6v as follows : 


6v = Tib^dx (2.6) 

Substituting equation (2.6) in equation (2.2) and simplifying, 
we get, 
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(2.7) 

2.2.2 D etermination of excess surface tension force F 

" ■" ' - ■'■■■"' ' - ' 1 ' '■ " ■ " s 

A curved surface is characterised hy two radii of cur- 
vature and R 2 which are given by the expressions, 

1/R^ = -y” /d+y'^)^^^ (2.8) 

1/^2 =i/[yd+y'^)^] • (2.9) 

The excess pressure (a R) at point P on the ellipsoid 

96 

(Fig. 2.2) is given by , 

aP = y(|^ + 1^) (2.10) 

where y is the surface tension of molten metallic particle. 

The component of force along the x-axis acting on a 
circular strip passing through P is therefore given by 

dF=APdAn.i.- (2.11) 

After substituting R^ and R 2 from equations (2.8) and (2.9) 
respectively in equation (2.10), the resulting aP and dA and 
n respectively from equations (2.3) and (2.5) are substituted 
in the equation (2.11). Then integrating the resulting expre- 
ssion from X = o to X = a, we get the excess force F along the 
x-axis given by 

Fg = PTcyb . 


( 2 . 12 ) 
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2.2,3 Determination of critical electric field 

Figure 2.3 represents the situation when two grains of 
identical radius r are stretched into ellipsoidal shapes under 
the application of an electric field. The net surface tension 
force acting in a direction opposite to electrostatic force is 

given by 

F = 2TtY(r-b) . (2.13) 

Xio 

At equilibriiom the electrostatic force F^^ and the net surface 
tension force (F ) must balance. Hence, 

•^Tiee^E^b^ = (2.14) 

Asstiming that the volume of the particle remains unchanged 
after stretching, it can be shown that 

b^ = r^/a (2.15) 


Prom figure 2.3 it can be shown that 


a=r+(sQ~s)/2 (2.16) 

Using equations (2.15) and (2.1 6) we get from equation (2.14) 
after simplification, 



The critical field (or switching field) is defined as the 
field for which the two ellipsoidal particles touch each other 
in equilibrium. So by putting s = o in expression (2.17), we 


get the critical field E^, given by 




Fig. 2.3 Two metallic grains stretched into 

ellipsoidal shapes by an electric field . 
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E = { 
m se_r 


^[{ 1 + - ( 1 + ]> 


( 2 . 18 ) 


2^2.4 Solution of one dimensional heat equation 


E^ has been derived on the assumption that the metallic 
particles are in the molten state. In our model the temperature 
rise of the metallic species is brought about by Joule heating. 
For delineating the corresponding temperature profile, within 
the sample we consider the following one dimensional heat 
equation ; 


^ ~ h 

dx 


where o be the electrical conductivity and 
conductivity and T is the temperature in 


( 2 . 19 ) 

K is the thermal 
degree Kelvin. 


The boundary conditions are : 


at X = o, 

at X = 1 

where 1 is the thickness of the sample and is the tempera- 

97 

ture at which the experiment is carried out . It has been 
45 

established for the silicate system under consideration that 
the electrical conduction takes place due to electron tunne- 
lling between the metallic particles. So, we have taken the 

47 

expression, for high field conductivity as given by 

oo 

cTjj = {qxjCD exp(-EyE) J dz z}* 

-V® 

Ec 

{ exp ( -z ) / [ 1 -exp ( — ■ g /g • ■ '£j') ] } 

o 


T = T 


exp 


tnd 


1=1 


exp 


( 2 . 20 ) 
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where ojj is the high field conductivity and k is the Boltsmaim 
constant. Ooc? and and are given by. 


<^00= ®^l'T'/[(sQ'+dQ’ (2 X s^* + 1 hj 


( 2 . 21 ) 


E^ = ojeis' + d^') 

O O O O 


=0 = 


( 2 . 22 ) 

(2.25) 


where 1* be the recombination length, y* 'tiie geometric 
factor, Sq ’ be the most probable separation between the 
two metallic grains of diameter d^’, e be the electric charge, x 
is a constant, h is the Planck’s constant, s be the separation 
between the two metallic granules and iP is the energy required | 
to generate a pair of fully dissociated positively and negativeHy 
charged metallic grains. 


For ^2^5"’^2^5 glass system the temperature dependent 


conductivity has been taken as 


44 


0 = cr exp(~9/kT) 


(2.24) 


where 9 is the activation energy. 

The heat equation (2.19) is nonlinear and no analytical solu- 
tion is available. Therefore, it has been solved numerically. 

2 . 2.5 Mumerical solution 

Equation (2.19) has been normalised, by putting y = x/l 
as follows ; 



33 


3 % = - o e2 

The houndhry conditions in this case become 


(2.25) 


® at y = o 

and 

* = ’^exp at y = 1 . 

OQ 

Fourth order Runge-Kutta-Nystrom method has been used to 
solve equation (2.25). Ihe normalised length has been divided 
into 10 equal parts with a step length h' =0.1. To use the 
R\mge-Kutta~Nystr6m method for second order differential equa- 
tion, it is necessarj’' to know the functional value and the 
first derivative at the same point. Since, in our case, we 
know only the functional values but not the first derivative, 
so to start with a trial value of slope at the boimdary has 
been assumed. In the case of silicate system, the conductivity 

expression (2.20) involves the integration and this has been 

/ ■ 99 

evaluated by using Simpson 1/3 rule . After computing the 
temperature distribution, we have compared the computed tempera- 
ture at h' =1 with the boundary value at h' = o. Depending 
on whether the computed temperature at h' = 1 is greater or 
less than the temperatirre at h = o, the trial value of slope 
is reduced or increased respectively. This process is continued 
until we get two slopes in subsequent iterations such that in 
one case the computed temperature at h' =1 is greater than the 
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boundary value and in the other case the temperature at h* = 1 
is less than the same. By using Regula~Palsi method, we have 
computed the temperature distribution, within an error of +0.1°K, 
The same technique has been used for a particular value of electric 
field to find the thickness within +0.1 i«n error such that the 
centre of the sample reaches the melting temperature of bismuth. 

The above scheme has been used to compute the temperature distri- 
bution for different combination of particle sizes, particle dia- 
meter to separation ratios, electric fields etc. In the computa- 
tion the appropriate values of k and s have been taken from 
references ([lOO], [44] and [45])- In. the numerical computation 
of the temperature profile, the effect of particle size on the 
melting point of bismuth to be described in the following section, 
has been incorporated,. The computer programme written in fortran 
language is given in Appendix 1 . i 

2.2.6 Dependence of melting point on particle size 

It has been observed experimentally^ ^ ^ that the melting te 
perature of bismuth particle decreases as the particle size de- 

j 

creases. The implication of this in our model is that the appliec 
electric field necessary to bring the meta3.1ic bismuth in the 
molten condition will be small for smaller particles. So the 
functional relation between the melting temperature and the 
particle size has been derived thermodynamically in the following ; 


section 



35 


2.2.6,. 1 Derivation of particle size dependent melting point 
expression 

Let a and p in figure 2.4 denote the solid and the 
liquid phases respectively under equilihrium at the melting 
temperature T of the solid phase. In such a case, it can be 
shown that under constant pressure. 


(S^ - S^)dT = -(dp^ ~ dp^), (2.26) 

where S denotes entropy, p the chemical photential, and the 
subscripts indicafe the phases. 

For a transfer of dn^ moles from a to p phase, the Gibbs 
free energy change at constant temperature and pressure is 
given by 

“It,P = ('V - ’ (2.27) 


where dA and dA^ represent the change of surface area in the 
two phases and y is the surface energy (we assume that solid- 
liquid and liquid -matrix interfacial energies have the same 
value ) . 

For the geometry assumed we can write 


= (2/r)\ 


(2.28) 


^-'^6 = 


HdR 


R dR - r dr 


■)Vpdnp, 


and 


(2.29) 
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where are the molar volxxmes of a and S phases respecti- 

<x p 

vely. At the start of melting, we can assume R equa- 

tion (2^29) then becomes 

¥e can also write 

4Ti;r^dr?p/?^ = 47iR^dR, (2.31) 


where terms on both sides represent the change of volume of 
the p phase as a result of the transfer of dn^ moles from the 
solid to the liquid phase. Assuming R 1. r, equation (2.31 ) 
simplifies to 


4R/dr =yv^ 


(2.32) 


Combining equation (2.30) and (2.32) we get 


2 

P a 


(2.33) 


Substituting equations (2.28) and (2.33) in eqiiation (2.27) and 


using the relation dn„ = -dn_, we get 

OC p 


d01j_p = (n^-lip)dn^ + % 




^ )dn 


\ 


a 


(2.34) 


At equilibrium, dG|g, p = 0. Therefore from equation (2.34) we 
can write 

P - 

diip - dn^ = - 2 ^ -■ ) dr . 


V - V 
p a 


(2.35) 
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For the transition from a to p phase, we can write 


AH - TCSp ~ S^) = 0 , (2.36) 

1 ♦ © f 

Sp - S„ = aH/T , 

where a H is the molar heat of fusion. Prom equatioiB (2. 26), 
(2.35) and (2.36) we get, 

^ dl = ^ (V - g . -~- ) dr. (2.37) 

■ r^ “ - V 

pa 

Integrating between the limits from (bulk melting point) 
to T and from r = octo r the following expression is obtained : 




In 


2 22L tf . —Ifi ) 1 

'J - A.V! r/ — 


m 


aH 


( 2 . 38 ) 






P ■ a 

The values for different parameters appearing in equation (2.38) 
have been taken from references( [102] and [l03]). 


2.2.7 Relation of computed fields to bulk field 

The computed critical fields E and E are assumed to be 

me 

■n T_ 

related to the fields and respectively calculated on the 
basis of thickness by the following expressions. 


E^ = (1 + d/sQ)E^ , (2.39a) 

E^ = (1 + d/sQ)E^ (2.39b) 
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where d and are average particle diameter and inter ~ 
particle separation. This is justified because of the large 


47 

glass^ metal » 


value of P T / P j. ratio 
2.3 Results and Discussions 


P being resistivity. 


Figure 2.5 is the plot of equation (2.38) showing the 
variation of the melting temperature of bismuth as a function 
of particle size. It is evident from the figure that the 
melting temperature decreases with the decrease in the radius 
of metallic particles. The implication is that for snaller 
particle sizes, the electric field necessaiy to reach the 
melting temperature of bismuth will be reduced. 

Figure 2.6 shows a typical temperature profile as obtained 
by solving one dimensional heat equation (2.19). Open circles 
in the figure represent the bismuth in silicate system with 
the radius of bismuth particle 250 = 1.4x10^ V/cm and 

s^/d = 0.0 1 . Closed circles represent bismuth in ^2^5”^ 2^5 
system with radius of bismuth particle 400 = 0.4x10^ V/cm 

and s^/d = O.Ol. The calculation is carried out at a tempera- 
ture around 160°C so that the theoretically computed critical 
fields can be compared with the experimental results. 

Figvires 2.7 and 2.8 represent the bulk critical memory^^^^^^^^^^^ 
switching field (equation 2,39a) as a function of different 
particle radii and separation to particle diameter ratios for • 
bismuth in the silicate system and ^2*^5'”^2^5 respectively. 



Temperature , K 


Fig 2 5 Var lation of melting temperature of bismu 
metal as a function of particle radius . 


Temperature 



Fig. 2.6 Temperature profile from solution of one 
dimensional heat equation in a sample of 
thickness 10 am . 



• • 
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It is evident from figure 2.7 that the threshold field in 

silicate system varies over two orders of magnitude 

(10-^ ~ 10^ V/cm) for a particular value of particle radius as 

the s^/d varies from 0.01 to 0.1. Similarly for vanadi^ml 

phosphate glass (figure 2.8) the threshold field varies over 

2 4 

two orders of magnitude (l0 - 10^ 7/cm) as s^/d varies from 

0,01 to 0.1. In table 2.1 values of critical field for 

"H 

melting the metallic bismuth and field needed for memory 
action are compared for different values of particle radii and 
separation to diameter (s^/d) ratios for a sample of thickness 
10 lan. It is observed that for parameters which can be realised 
in practice, the E° is in general greater than e“. It is, 

HI W 

therefore, concluded that the threshold field for the ^sterns 
under investigation will essentially be controlled by E^. 

In table 2.2 the threshold field values predicted by the 

I 

present model are compared with the experimental results 
reported to-date for some silicate and vanadium phosphate I 

I 

glasses containing metallic bismuth. It is seen that the i 

i 

computed critical electric fields are of the same order of | 

magnitude as observed experimentally. Ihe agreement between 
these values are reasonable considering that we have made a 
number of simplifying assximptions in our analysis. 



'x‘rble 2.1 Critical fields for melting (E^) and memory action (E^) in gl 
matrix containing bismuth granules 
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Table 2.2 Comparison of theoretical and experimental 
switching fields 


Glass Composition 
(mole ia ) 

Maximum 

Diameter 

i 

E^/Exp 

.m 

( ¥/ cm ) 

EVcalc 

m 

(V/cm) 

643102, 18 B20^ 

500 

1 ,2x10^ 

1 . 1xl0^( Sj^/d=0 . 03 ) 

8 Bi20^,l0 lTa20 



2,4xl0^(syd=0.05) 

3.9xl0'^(sydd0.07) 

5.6xl0^(sQ/d=D.09) 

807205,1 5P2O5 

800 

4.0x10^ 

1 .0xl0^(sQ/d=O.0l ) 

5 Bi205 



2.8xlO^(sQ/d=0.02) 

5.1xl0^(s^/d=0.03) 1 
7.8xl0^(syd=0.04) 1 




CHAPTER 3 


CURRENT -^-YOITAGE CHARACTERISTICS OF VAHAIHIM PHOSPHATE 
GLASSES CONTAIHIFG BISMUTH GRAHULES 

3 . 1 Introduction 

In Chapter 2, the threshold field for memory switching has 
been calculated in the case of glasses containing microgranules 
of bismuth dispersed in them. The I~V characteristics of the 
present systems of glass~metal composites show negative resis-” 
tance region'^^’ . In the present chapter, the particle 
stretching model which has been developed in Chapter 2 has been 
used to compute the I-T characteristics of vanadium phosphate 
glasses containing bismuth, especially in the negative resi- 
stance region prior to the on-state as a function of particle 
diameters and particle separation to diameter ratios. The 
effect of specimen thickness on switching behaviour has been 
delineated both theoretically and experimentally. Also, the 
effect of temperature on the resistance of the samples in the 
memory state has been studied. 

3-2 Theoretical Consideration 

3.2.1 Calculation of I-V characteristics in the negative 
resistance region 

The threshold field as calculated in Chapter 2 for memory 
action corresponds to the beginning of the negative resistance 
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region in the I-V plot of the glass-cietal composite. According 
to the present model, initially within a portion of the sample 
at the centre where the temperature attains the value corres- 
ponding to the melting point of bismuth, elongated metal grains 
touch each other and thereby form a highly conducting region. 
Due to higher current flow more heat is generated and the 
conductive region grows from centre outwards. Such a regenera- 
tive process results in negative resistance characteristics and 
ultimately the system atta,ins the highly conducting on-state. 


For calculating the current -voltage (I-Y)values in the 
negative resistance region we consider in figure 5.1 the' 
situation within the specimen at some instant of time after 
the application of the critical switching field. In this 
figure 1^ indicates the width of the conducting region 
(referred to as metallic region in subsequent discussion) where 
molten granular metal particles are interconnected and 1^ 
represents the overall width of the remaining portion (referred 
to as dielectric region in subsequent discussion) of the sample. 


The resistance of the metallic region is given by, 

central DfBRAR" 

1 L i. 7 k 




( 3 . 1 ) 


where is the resistivity of the specimen in the on-state. 

Experimentally observed in the on-state has been used in the 


present analysis. A is the electrode area of the specimen, 




^ > 

*' m I 

— i - 


k- 


Fig. 3.1 Schematic representation of a glass - 
metal composite after the application 
of a critical field . 

Irr, ; width of metallic region 
: width of dielectric region 
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Since resistivity of tiie glass-metal particulate system 
is dependent on temperature, to find the resistance of the 
dielectric region one has to take into account the temperature 
profile within the sample after the application of the critical 
electric field. The dielectric region has been divided into 
n equal parts each of width w so that the resistance of this 
region is given hy 



n 

S 

i=1 



(3.2) 


where denotes the resistivity of ith strip having a mean 
temperature T^ and 

w = .1^/n (3.3) 

In our calculation we have chosen n suitably such that at 
every step of metallic region formation, the value of w is kept 
constant. The temperature distribution in the sample has been 
obtained by solving the time -dependent heat equation to be des- 
cribed in the following section. The resistivity has been 
evaluated from the equation (2.24) in Chapter 2. The pre- 
exponential factor and the activation energy have been obtained 
from experimental data of the present vanadium phosphate system 
containing bismuth^^. 

If the ciirrent passing throiigh the circiiit consisting of 
the sample and the series resistance R^ ohms be denoted by 1^ 
amperes, we can write. 
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and 




E._ 


m 


+ Rj + R 
d s 


(3.4) 


'5.5) 

where is the applied voltage across the circuit. 

For delineating the temperature profile within the sample, 

the rate of heat generation per unit volume needs to be evalua- 

2 

ted. The latter is given by o^E , where is the conductivity 
of the specimen and E the electric field. However, after the 
onset of switching, in the present model, the sample consists 
of two regions viz. metallic and dielectric. The average 
values of field and conductivity to be used for computing the 
rate of heat generation per unit volume are therefore obtained 
as follows : 


We can assume that after the coalescence of neighbouring 
metal granules near the specimen centre, the potential drop 
in the dielectric region occurs only in the glass matrix inter- 


vening the particles 
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With this assumption, the effective 


field is then given by 


E 


eff 




(3.6) 


where, s^ is the average interparticle separation and d the 


average diameter of metal particles. 


Substituting and R^ respectively from equations {3.4) and 
(3.2) in equation (3.6) we get. 
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E’ 


eff 


v^ppd + (d/s^)) 


n 

W V 

I i=1 


T 


n 


(3.7) 


(R + R + r^ Pm ) 

^ m s A • > T. ^ d 
x=1 1 


Similarly, the effective field, in the metallic region is 


given hy. 




I R 
3 m 

^m 


(3.8) 


Prom equations (3-1 ), (3.4) and (3.8) we can write. 


pin 

eff 


7 P' 
-.aPP., ..-S- 


+ Ss + i 


n 


(3.9) 


The average electric field in the glassHnetal particulate 

3, V 

system at any time after the onset of switching is then given 

( 3 . 10 ) 


B 


av 


= (1, + 1 #„.)/i 

^ d eff m eff'^ 


The average conductivity cr is given hy, 

3 V 


air 


^ f ^ + J^) L 
d m 


(3.11) 


where a. and a refer to the conductivities of the dielectric 
dm 

and the metallic region respectively. 

The present analysis has been carried out with the applied 
field which correspondsto the critical field for memory action 
^m the expression (2.18) in Chapter 2- is related 

to the bulk field E^, calculated on the basis of sample thickness 

by. 
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= (1 + <i/s^)E^ (3.12) 


where 

Prom equations (3.12), (3.13) 
Chapter (2). ¥e get 


(3.13) 

and equation (2.18) of 


\pp : ^1+d/s, ^ ^ 




(3.14) 


o o 

In the above expression, e,e^,Y,r have the same meaning as 

given in Chapter 2. To find the expression for 7 the 

app 

voltage applied to the circuit, we have neglected the poten- 
tial drop across the resistance in series with the sample. 

This is because of the fact that the magnitude of series 
resistance is orders of magnitude less than that of the sample 
in the off-state. 


3.2.2 Numerical scheme for I-V computation 

The normalised thickness is divided into ten equal parts 
and the grid points are denoted by 1,2, 11 respectively. 

First of all the resistance of the sample is calculated from 
equation (3.2) corresponding to ambient temperature. From 
equations (3.4), (3.5) and (3.14) the potential drop and the 
current passing through the sample are then evaluated. These 
correspond to the starting point of the negative resistance 
region. In our calculation, we have neglected the temperature 
rise, due to Joule heating in the pre -switching region. The 
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time -dependent heat equation is then solved numerically to find 

a temperature distribution such that the sample centre (grid 

point 6) has attained the melting temperature of bismuth 

corresponding to the average particle size (equation (2.58]t 

Chapter 2). Solution is then carried out corresponding to the 

time when the grid points 5 and 7 reach the melting temperature 

of bismuth. The molten particles in the region lying between 

the grid points 7 and 5 will be stretched forming the metallic 

region, because the electric field has the critical value. 

Va,lues of and R^ are then calculated corresponding to this 

state from equations (5.1) and (3.2) respectively. Substituting 

these in equations (3.4) and (3.5) the I and Y values are 

s s 

calculated. The values of E and a are then estimated from 

av av 

equations (3.10) and (3-1 l) respectively, which are then used 

in the next step of the computation. This procedure is 

repeated and I ,Y values are found at every step till we 
s s 

come to the one which involves the pair of grid points 2 and 
10 respectively. The latter give the last but one point in 
the negative resistance region. The final point corresponds 
to the situation when the entire sample forms the metallic 
region thereby attaining the on-state. 

3.2.3 Relaxation of temperature 

To explain the threshold switching behaviour of the thick 
sample, as we shall see later, it is necessary to calculate the 
time required for the temperature at the centre of the sample 
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to go telow the melting point of bismuth after the applied 
electric field is reduced to zero. The relaxation time (t) 
is defined as the time which elapses after the removal of the 
electric field for the on-state temperature at the centre of 
the sample to attain the melting temperature of bismuth. The 
decay of temperature distribution in a rectangular slab after 
normalisation with respect to thickness is given by, 


df K d~'T 


(3.15) 


where k , P' and c respectively denote thermal conductivity, 
density and specific hea-t respectively. 


The boundary conditions are : 


T(o,t) = T(l,t) ^ 0 
T(x,o) = f(x) for o <x < 1 


( 3 . 16 ) 


where f(x) is the temperature distribution within the 
sample when it attains the on-state. 

98 

Solution of equation (3.15) is given by 

oc - xjt 

T(x,t) = Z [ B e sin(n'n;x)] (3.17) 

n=1 ^ 

where , = ~ ( ^ (3 .18) 

1 

®n ~ 2 J f(x) sin(n'n;x)dx 
0 


and 


(3.19) 



59 


The solution of equation (3.17) has been worked Out numerically, 
since only numerical values of temperature distribution function 
f(x) are known in the on-state. 


3.2.4 numerical solution 

3 . 2 . 4 . 1 numerical solution of time dependent heat equation 

The time dependent heat equation which is solved to find 
the temperature distribution within the sample for computing 
the I“V characteristics in the negative resistance region is 
given by 

p.c .^ 3 , ic^+ ar{T)E^ . (3.20) 

dt* dj 

The boundary conditions are : 


T(y) = °< y # at t* = o 


T(o) = T(l) = T^^p 


at t* > o 


For ease of numerical computation, the equation (3.20) has 
been normalised with respect to thickness and for some time t^ 
by setting x = y/l and t = t*/t^ and the normalised equation 
becomes. 


P c o>T 


JL iZl 

i 2 .^2 

1 ffX 


+ o(T)E' 


The boundary conditions become , 


(3.2T) 


T(x) = Tg^p for 0 <x <1 at 


T(o) = T(l) = T 


exp 


t ■ o"" 

at t > 0 
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For solving equation (3.21 ) we have used Crank-ITicolson 
implicit finite difference formula'^^^. The advantage of 
implicit finite difference scheme is that it gives unrestri- 
cted stability for p/h^ ratio^*^^, where p and h are respecti- 
vely the step lengtte in time and x direction. So we have 
chosen a step length p =0.1 and h = 0.1 in our computation. 

In Crank-Nicolson implicit formula the time derivative is 
replaced by forward difference and the second order derivative 
in the x direction is replaced by a weighted sum of second order 
central difference at t and t+p. Applying these to equation 
(3.21 ), we get after simplification the following finite 

<• 

difference formula at a particular time t, 

.. Jl_ T + fiLa-lc . _ JL- T 

x-h,t+p ^ t^p ^ 2 ^ x,t+p x+h,t+p 

2 

— r — m 4, f sL IT 

- ‘ t,p 

+ i^g(- ^ -*- it£ (3.22) 

In equation (3.22) the conductivity has been evaluated at a ■ 
mean temperature of and Writing equation (3.22) 

for all the grid points along the x -direction we obtain a (9x9) 
tridiagonal matrix. For finding the temperature distribution 
from the above equations the Gauss matrix iteration^^ has been 
used. It can be shown from the above ( 9x9 )tridiagonal matrix 
that each diagonal element is greater than the absolute sum of 
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the corresponding off-diagonal elements. This satisfies the 
condition for convergence of Gauss iteration scheme. 


The following steps are involved to calculate the tempera- 
ture distribution at different time intervals and to compute 
the time necessary to reach the melting temperature of bismuth. 


Step 1 

Step 2 
Step 3 

Step 4 

Step 5 

Step 6 

Step 7 

Step 8 


A temperature distribution is assumed to start the 
Gauss iteration, t = o and p = 0.1 are set, and a 
value of t^ is assiomed. 
cr is evaluated corresponding to time t. 

Temperature distribution is computed using cr as 
evaluated in step 2. 

cr is evaluated at the mean value of temperatures 
evaluated at time t and the temperatures evaluated 
in the last iteration. 

Temperat\xre distribution is obtained using cr as 
evaluated in step 4. 

Steps 4 and 5 are repeated until the difference of 
temperature distributions between two successive 
iterations is less than 0.1*^K. This gives the tem- 
perature distribution at time t+p, 
t = t+p is set and steps 2 to 6 are repeated till 
t = 0.9. 

The difference between the temperature at the centre 
of the sample and the melting point of bismuth is 
computed. 
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Step 9 : 

Step 10; 

Step 1 1 : 

Step 12; 

Step 13; 


Depending on whether the difference is positive or 
negative t^ is decreased or increased. 

The value of t^ in step 1 is replaced by the new 
value of t^ of step 9. 

Steps 1 to 10 are repeated until we get two values 
of t^ such that in one ca.se the difference as evalua- 
ted in step 8 is positive and in the other case it is 
negative . 

Using method of bisection new value of t^ is evaluated 
and t^ in step 1 is replaced by this new value of t^ . 
Steps 1 to 8 and then step 12 are repeated until the 
difference as computed in step 8 is within + 1°E and 
the time which is evaluated in step 12 is the reqiiired 
time . 


The programme for the above nmerical computation is given 
in Appendix 2. 

3-2. 4-2 Numerical solution of temperature relaxation 

The coefficient (equation (3.19))haS’ been 

evaluated by using Simpson l/3 rule using the numerical values 
of f(x), the temperature distribution as obtained from the 
solution of time dependent heat equation for different particle 
diameters and separation to diameter ratios. Using the method 
of bisection, the relaxation time as defined in Section 5.2.3 
has been evaluated numerically for thin as well as thick samples. 




63 


3 .3 F.X -per i mental 

3.3.1 Preparation of glass 

A glass latch of composition (mole 5i5 ) 

has been prepared from reagent grade chemicals. ^ 2^5 added as 
(NH ^)2 HPO^ and "^ 2*^5 ®^ 2^3 introduced as their oxides. 

The mixture is melted in a pure crucible at temperatxire of 

700°C for half an hour. Then the molten glass is cast in an 
alxmiinium mold and annealed at 100°C for half an hour and left 
for furnace cooling. 

3-3.2 Sample preparation 

3 . 3 . 2 . 1 Preparation of thin film 

44 

The thin film of glass has been prepared as follows : 

The liquid droplet of molten glass is collected by dipping a 
stainless steel tube in the crucible containing the molten 
glass. The droplet has been blown into a large bubble by 
passing an inert gas e.g. nitrogen/ argon through the tube from 
the other end. The thin film samples are collected from the 
resulting glass bubble. 

3. 3. 2. 2 Preparation of thick film 

The bulk glass has been thinned by successively grinding 
and polishing in SiC powder. Since vanadium phosphate glaisses 
are readily attacked by water, the grinding and polishing are 
carried out in the presence of acetone. Samples of thickness 
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a few hundred microns are prepared hy this method and are 
termed thick films. 

3.3.3 Differential thermal analysis (DTA) 

The bulk sample is ground and then sieved out in the 
range of ~1 4 mesh to 20 mesh. This particle size is large 
enough so that the effect of surface crystallisation can be 
neglected in DTA . To the ground sample is added aroxmd 
30 io of calcined ■^2^3 avoid sintering of the glass par- 
ticles. The DTA has been carried out in a MOM Derivatograph 
Hungary, at a heating rate of 7°C/minute. 

3.3.4 X-ray analysis 

Both heattreated and virgin glasses have been used for 
X~ray studies- The heat treatment of sample has been carried 
out in a wire wound furnace at 359°C for a period of one hour. 

The glass sample is ground and sieved to get powder of size 
below 400 mesh- The X-ray diffraction pattern of the powdered 
sample has been taken using General Electric XRD 5 Diffractometer 
The monochromatic radiation used is Cu The interplanar 

spacings have then been calculated from the intensity vs 
2d plot of the diffraction pattern by using the relation, 

2d, Sind = n A (3.23) 

hkl 

where and ©■ are the interplanar spacing. X-ray wave- 

length and the diffraction angle respectively. 
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3.3.5 Electron microscopic analysis 

Electron microscopic analysis has been carried out for 
both heat treated and virgin glass samples. The electron 
microscopic sample has been prepared as follows ; Eirst, 
the sample is thoroughly ground in an agate mortar and then 
dispersed in acetone, so that the very fine particles, needed 
for transmission electron microscope remain in suspension. A 
few drops of acetone containing finely suspended glass particles 
is poured on the carbon coated electron microscope (EM) grid. 
One drop of collodion is then poured on deionised distilled 
water to form a very thin film of collodion. The EM grid 
containing the powdered glass sample is then put on the film 
so that the sample particles face the collodion film, lastly, 
the EM grids are fished out by clean glass slides. The 
sample is, therefore, sandwitched between the carbon and 
callodion thin films. The electron micro structure and the 
selected area diffraction of the samples at lOO kV have been 
taken in a Phillips EM 301 Electron Microscope. 

The diameters of the diffraction rings have been measured 
from the negative by using a graduated eye-piece with an 
accuracy of 0.1 mm. The interplanar spacings are calculated 
from the relation given below ; 

hkl 


Camera constant = r d. 


(3.24) 



the 


where r is the radius of the diffraction ring and d ^ ^^ 
corresponding interplanar spacing. The camera constant has 
been determined from the standard gold diffraction pattern 
obtained at 100 kY. 

3.3.6 Electrical measurement 

The I-”Y characteristics of both thick and thin samples 
have been measured. The resistance of the thin film sample in 
the memory state has been measured as a function of temperature. 
To carry out the electrical measurement, the samples have been 
coated with gold to form electrodes. This has been done by 
evaporating pure gold in vacuum, 

3.3.6. 1 Electrical circuit 

Figure 3.2 represents schematically the electrical circuit 

used for measurements of I-Y characteristics of different 

samples. The sample and the standard resistance are in series 

with a 30 Yolts 2 amperes DC regulated power supply (Aplab, 

Model 7152). The potential drops across the sample and the 

standard resistance are fed to the X-axis and Y-axis of the 

X-Y recorder respectively (No. 2000 made by Digital Electronics 

Ltd., India). Hence, the plot on the X-Y recorder corresponds 

to the I -Y characteristics of the sample. The same circuit 

has been used to measure the resistance of thin samples in the 

.0 

memory state as a function of temperature from 23 C to about 




Fig. 3.2 


Schematic representation of electrical 
circuit . 
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107 

5-5. 6, 2 MeasuriBff cell 

Figure 5.5 shows the schematic view of the cell used 
for electrical measurements. This consists of a double 
walled pyrex glass Jacket J, 500 mm long and having a diameter 
of 50 mm. and are two gold plated copper electrodes 

being placed co-axially with the metallic hollow cylinder S. 
They are separated from the metallic hollow cylinder S by two 
teflon discs T. The lower electrode is fixed while the 
upper electrode can move up and down. The electrodes are 
connected to the BEC connectors B with shielded wires. The 
sample can be placed between the two electrodes by inserting it 
through a. rectangular hole 0 in s. On the other side, a copper 
constantant thermocouple is inserted throvigh a small hole H. 
The sample holder assembly is suspended with the help of three 
brass rods R brazed to the upper copper flange w f with a 
rubber gasket G in between, H.C. is a heating coil wound on a 
porcelain tube. This is connected to a 50 volts 2 amperes 
regulated DC power supply through a proportional temperature 
controller. The copper constantant thermocouple TO serves 
the purpose of measuring temperature as well as sensing ele- 
ment for temperature control. By this set up the temperature 
can be varied from room temperature to about 1 50°C with an 
accuracy of temperature control + t^C. The inside of the 
assembly can be flushed with dry inert gas. To make the inside 
of the assembly moisture free silica gel dessicant D is used. 
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5.4 Results and Discussion 

3.4.1 Switching characteristics of thin samples 

Figures 3.4(a) and (b) show the theoretical chara- 
cteristics obtained for a film of thickness 10 and electrode 
area 6.32:10"’ m for different particle diameters and particle 
separation to diameter ratios. In figure 3.5 are shown the 
theoretical temperature distribution in a sample at different 
instants of time after the application of the critical field. 

It is seen that the time taken for the sample centre to attain 
the melting temperature of bismuth corresponding to a, particle 
diameter 700 2. and separation to diameter ratio 0.03 is 
1.72:10 second . In column 3 of table 3.1 we give the calcu- 
lated time tj^ for the entire sample to reach the melting tem- 
perature of bismuth for different values of d and s^/d ratios- 
fhe kinetics of switching at room temperature is controlled by 
the heating step. However, we will show in Chapter 4 that the 
stretching kinetics is much faster. 

In figure 3.6 the I-V characteristics obtained experimen- 
tally on a sample of identical geometry to that considered above 
are shown. It is noted that the sample switched to the memory 
state when potential drop across the sample is aro^md 6,3 Volts. 
Also plotted in the figure are some theoretical I ~V curves taken 
from figure 3.4(a) for particle diameter values of 300 % and 
400 2 . and s^/d ratios equal to 0.02 and 0.03, The agreement 
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Fig. 3.4(a) Theoretical I-V characteristics for a glass- 
^ ^ ^ ^ ^ 80V2O5 15P2O5 

5 Bi 203 (mole %) for various particle diameters 
Sample thickness = 10 |U.m 

area = 6.3 x 10 m^ 
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Theoretical I-V characteristics of glass - 
.metal system of composition 8 OV 2 O 5 I 5 P 2 O 5 
5 Bi 203 (mole %) for various particle diameters 
Sample thickness = 10 fim 
Electrode area = 6.3 x 10 m^ 
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Table 3.1 Calculated ps,rame-ters for switching in a 
glass of composition 80720^!> 15I’2*^5>5Bi2*^3 
(mole ^ ) for different d and s^/d values 


Separation 
to diameter 
ratio 

So/d 

Switching time t^^ 
controlled by- 
heating step 

xlO'^ sec. 

1=1 O^m 1=45 lan 

Relaxation time x for 
temperature profile 

xio"^ sec. 

I=l0iim l=45iffii 

0.02 

1.32 

0.69 

0.11 

0.42 

0.03 

0.72 

0.45 

0.10 

0.10 

0.04 

0.50 

0.34 

0.10 

0.10 

0.05 

0.36 

0.27 

0.09 

0.10 

0.06 

0.28 

0.22 

0.09 

0.10 

0.07 

<M 

• 

O 

0.20 

0.08 

0.10 

0.02 

7.39 

5.73 

0.15 

2.27 

0.03 

3.88 

3.09 

0.15 

1 .83 

0.04 

1 .82 

2.00 

0.15 

1.51 

0.05 

1.36 

1.40 

0.17 

1 .32 

0.06 

1 .07 

1.06 

0.16 

1.14 

0.07 

0.89 

0.85 

0.16 

1.03 

0.02 


11.88 


2.50 

0.03 

5.01 

6.90 

0.15 

2.32 

0.04 

CO 

00 

(M 

4.73 

0.16 

2.11 

0.05 

2.05 

3.45 

0.16 

1.97 

0.06 

1.59 

2.65 

0.16 

1-79 

0.07 

1.29 

2.05 

0.16 

1.59 


contd, . . 
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0.02 


21.72 

- 

2.55 


0.03 

8.49 

11.77 

0.14 

2.44 


, 0.04 

4.15 

8.20 

0.14 

2.35 

600 

0.05 

3.41 

6 .08 

0.14 

2.22 


0.06 

2.12 

4.76 

0.15 

2.11 


0.07 

1 .72 

3.80 

0.16 

2.02 


0.02 

- 

56.67 


2.55 


0.03 

18.46 

19.50 

0.14 

2.50 


. 0.04 

5.86 

12.55 

0.14 

2.41 

700 

0,05 

3.69 

9.20 

0.14 

2.32 


0.06 

2.72 

7.33 

0.15 

2.27 


0.07 

2.16 

5.91 

0.16 

2.19 


0.02 

— 

~ 

— 



0.03 

- 

- 

- 

- 


0.04 

8.60 

19.10 

- 

2.47 

800 

0.05 

4.82 

13.27 

0.14 

2.41 


0.06 

3.42 

10.25 

0.14 

2.35 


0.07 

2.65 

8.28 

0.14 

2.11 


0.02 


— 




0.03 

- 

- 


- 


0.04 


30.10 


2.50 

900 

■ 0.05 

6.36 

18.88 

0.14 

2.44 


0.06 

4.27 

13.95 

0.14 

2.38 


0.07 

3.23 

11.10 

0.14 ' 

2.35 


0.02 

— 

— 




0.03 


— 




0.04 

- 

67.95 

, . 

2.53 

1000 

0.05 

8.63 

29 . 1 1 

0.14 

2.50 


0,06 

5.31 

19.58 

0.14 

2.44 


0.07 

3.90 

14.90 

0.14 

2.41 
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Fig. 3.6 Comparison of theoretical and experimental I“V curves 
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between the two is found to be reasonable. Pigixres 3.7(a) and 
(b) gi-ye the electron micrograph and the selected area electron 
diffraction pattern of the sample respectively. The particles 
are observed to have diameters varying from 50 S. to 450 
In table 3.2 the interplanar spacings calculated from 

the diffraction pattern are comparod with the standard 

1 0R 

values for bismuth metal. The agreement between the two 
sets of values indicates that the particles comprise of metallic 
bismuth. 

3.4.2 Switching characteristics of thick samples 

The effect of specimen thickness on the switching behaviour 
has been investigated by measuring the I~V characteristics of a 
sample of thickness 356 Ho switching has been observed at 

room temperature and upto an applied voltage of 60 Volts. When 
the sample is however heattreated at 360*^0 for an hour it shows 
threshold switching with a threshold voltage of around 26 Volts. 

Also the resistance of the heattreated sample is about an 
order of magnitude lower than that of the untreated one. The 
heattreatment is believed to cause further pre-cipitation of 
metallic bismuth in the glass matrix thereby bringing about 
coalescence of larger particles which effectively gives rise to 
the formation of metallic region in a major portion of the glass 
metal particulate system. This is ' substantiated as follows ; 
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Table 3.2 Electron diffraction data for a 
glassHnetal composite 


'^hkl calculated 


from experimental data 

i 


standard d ^^ .^ 
bismuth, metal 


values of 
108 


i 


3.30 

3.28 

2.89 


2.18 

2.27 

2.07 

2.03 

1 .86 

1 .87 

1 .64 

1.64 

1.30 

1.31 

1.11 

1.12 
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Figure 3.8 shows the DTA curve obtained for the as-cast 
glass. The exothermic peak at 360°C indicates that some 
crystalline phase is precipitated within the glass matrix at 
this temperature. In figure 3.9 are shown the X-ray diffra- 
ction patterns of the as-cast glass and of the sample heated 
at 359°C for an hour respectively in the range of 26- values 
where the most intense peak of the heattreated sample occurs. 
Figures 3.10(a) and (h) give the electron micrograph and 
the selected area diffraction pattern respectively of a 
heattreated sample. Table 3.3 summarises the X-ray and 
electron diffraction data of the same. It is evident from this 
table that the precipitated phase consists of metallic bismuth. 
Also from figure 3.9, it is obvious that the amount of metallic 
bismuth present in the heattreated sample is higher than that 
in the untreated sample. From figxire 3.10(a) we conclude that 
the bismuth particles in this sample have diameters in the 
range 50^ to 700 2.. 

Figure 3.11 gives the theoretical 1-7 curves obtained for 

2 

a sample of thickness 45 1->m and electrode area iSxlO m for 
different particle diameters and separation to diameter ratios. 
The above thickness value has been chosen such that the experi- 
mental data match closely the theoretical results, for particle 
diameter values 400 2 . and 500 2 . and s^/d value equals to 0.03» 
as shown in figure 3.12. From the discussion given above the 
effective thickness should be an order of magnitude smaller than 
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Table 3.3 X~ray and electron diffraction data for 
a glass-metal composite heated at 359®C 
for 1 hour 


'^khi 

calculated from 

electron 

diffraction 

2 


calculated from 
X-ray diffraction 


Standard d,, t values 
hkl 

of bismuth metal 


i 


i 


— 

4.39 

- 

- 

3.37 

3.28 

- 

2.90 

- 

2.56 

- 

2.39 

2.23 

2.18 

2.27 

1 .84 

- 

1 .87 

- 

1 .50 

1.49 

1.42 

- 

1.44 

1 .29 

— 

1 .28 
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ig. 3.12 Comparison of theoretical and experimental I- V 
curves . 

Sample thickness = 45 jxm 
Electrode area = 16 x 10 ^m^ 
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the actual thickness (556 ^aa in our case). Our estimated 
value (45 therefore seems reasonable. In column 4 of 
table 3.1 we give the calculated time required for the entire 
sample to reach the melting temperature of bismuth for various 
d and s^/d values. The time t^^ for the thick sample is 
found to be in general greater than that for the thin sample 
which is expected from physical reasons. 

In contrast to thin sample the threshold switching 
behaviour of thick sample can be explained from the calculated 
values of relaxation time t which has been defined earlier 
(Section 3.2.3) as the time required for the temper atxrre at 
the centre of the sample to go below the melting temperature 
of bismuth after the applied electric field is reduced to 
zero. The relaxation time has been computed by the procedure 
outlined in Section 3.2.3 and Section 3. 2, 4. 2 for different 
d and s^/d values and for both the sample thicknesses under 
investigation. Fig'ure 3,13 shows a typical calculated tem- 
perature profile as a function of time for a sample of thick- 
ness 10 ima with d and s^/d values of 400 % and 0.03 respecti- 
vely, The results are summarised in the last two colimins of 
table 3,1. It is evident from these data that -x values obtained 
for the sample with a thickness of 45 are an order of magni- 
tude higher than those of the thinner one (10 jaa). If the 
rate of removal of the electric field is faster as compared to 
the relaxation time r, it is likely that the stretched particles 



Temperature , 



Fig. 3.13 Theoretical relaxotion of temperature profile 

within 0 sample of composition 80V205 15P205 
5Bi203 (mole %) after the removal of the 
electric field. Sample thickness = 10 u,m . 
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will regain their original shape due to the surface tension 
forces still effective in the molten state of metal granules. 
This will break the interconnectiors between the conducting 
particles and the sample will go back to the off -state. The 
occurrence of threshold switching in the case of the thick 
sample investigated here is therefore believed to arise due to 
a higher value of % pertaining to the specimen. 

3.4.3 Temperature dependent memory state resistance 
variation 

Figure 5.14 shows a typical variation of resistance in the 
memory state of the blown thin film as a function of temperature 
over the range 20^0 to 110°C. Before switching the off -state 
resistance of the sample is 1070 Ohms. It goes to the memory 
state via a negative resistance region at an applied voltage 
of 9.5V at 23°C. The resistance in the memory state is 8.0 Ohms 
at the same temperature. It is found from figure 3-14 that the 
on-state resistance remains constant with temperature. At 
107°C the sample goes to the off-state of resistance about 
8000 Ohms with an application of 0.5 volt . This behaviour can 
be explained as follows. According to particle stretching 
model, if the particles are in molten condition in the on-state, 
the removal of electric field causes the stretched particles to 
go back to their original shapes due to surface tension. Thus, 
the interconnections of the metallic chains are broken and a 
high resistance off -state is achieved. This condition is 
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realised by heating the sample in the memory state to 107°C 
along with an applied voltage of 0.5 Volt. It is believed that 
the role of the applied voltage is to help achieving the melting 
temperatiire of bismuth. 

A similar behaviour has been reported in the case of 

silicate glass system containing bismuth at high temperature 
51 

alone . 

When the interconnected chains in the memory state between 

the molten bismuth particles are broken due to the effect of 

surface tension force, the resistance of the sample increases 

and at the same time the increased resistance of the sample 

causes the ciirrent to decrease. Thus more the breaking of 

the interconnected chains the greater is the resistance of the 

sample and therefore the potential drop across the sample 

increases along with a reduction of current giving rise to a 

voltage controlled nega.tive resistance region in the I-V 

characteristics. Figure 3.15 is the I-V plot of another sample 

in the on-state at 78°C. The sample goes to off-state at 

0.18 Volt through a voltage controlled negative resistance 

region. Similar voltage controlled negative resistance region 

has been observed in the insulating thin film sandwitched 

between two metallic electrodes due to thermal rupture of 

109 110 

conducting filaments ’ . 

A n-umber of samples has been studied which show qualita- 
tiyely same behaviour as described above. 



Current (amp xIO 



Potential drop (volt) 


Fig. 3.15 I-V plot of 80V2O5-15P2O5-5Bi2O3 (mole %) 
gloss switching from on to the off -state 
at 78°C. 



CHIPTER 4 


KINETICS OP SWITCHIN& IN VANADIUM PHOSPHATE GLASSES 
CONTAINING BISMUTH GRAITOLES 


4 . 1 Introduction 

In Chapter 2, we have developed a model. This predicts 
the critical switching fields which are of the same order as 
those obtained experimentally. In Chapter 3 with the help 
of the model developed in Chapter 2, the I~V characteristics 
especially in the negative resistance region have been computed. 
There it has been assiimed that once the temperature reaches the 
size dependent melting point of bismuth interconnected chains 
are formed because of the presence of critical electric field. 
This implies that the time required for two physically proxi~ 
mate grains to touch each other in molten condition is much 
less compared to the time required for the temperature profile 
to reach the melting temperature of bismuth. In this chapter, 
we present the solution for the dynamics of stretching and 
compute the time required for two physically proximate grains 
to touch each other in molten condition after the application 
of critical switching field. 

4.2 Theoretical Consideration 

When an electric field is applied to a glass-metal 
composite, the metallic particles in the molten state will be 
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stretciied as described in Chapter 2. Figure 2.5 of Chapter 2 
shows, at some instant of time, the state of two proximate 
grains during stretching into ellipsoidal forms from spherical 
shapes. The problem is to find the time required by each 
sphere to elongate by a distance 8 ^/ 2 , so that the particles 
can join up. It is noted that the boundary of the ellipsoid 
changes during stretching. This moving boundary fluid flow 
problem has been analysed with the help of well known fluid 
flow equations. 

The following assumptions have been made in our analysis 
for simplicity of calculations. 

1) The molten liquid is inviscid i.e., the viscous drag 
has been neglected compared to surface tension. 

2) The flow is irrotational and the liquid is incompressible. 

5) The flow problem has been solved in two dimensions. 

4,2.1 Determination of potential 9 ^ 

To determine the force acting on the deformed granule at 
a distance x from the origin, we use the same method as des~ 
cribed in Section 2 . 2.1 of Chapter 2. By displacing the 
ellipsoid by 6 x (figure 4.1) against the force F , in the 

k 

negative direction of x, we can write, 

Fg^l 6 x = i ee^ B^ 6 ¥ (4.1) 

where 6V, is the volume as cross-hatched in figure 4.1. 



Y| E 

^ 



F'lg. 4.1 Configuration of a half ellipsoid 
stretched under the action of an 
electric field , 



97 

Its value is obtained by integrating the equation (2.4) of 


Chapter 2 from x = o to x = x and we get, 


^2 5 

6? = -rtC— 6x 
a 


(4-2) 


where a and b are semi-major and semi -minor axes of the 
ellipse respectively. 

Substituting equation (4.2) into equation (4.1) and after 
simplification, we get. 


^el I ^ 


ee^ e2(^)x2 
a 


(4.3) 


This force will act on the mass of the ellipsoid between 
X = o and x =x. The volimae (V) of ellipsoid between x = o 
and X = X is given by. 


7=1 It y'^ dx 
0 


(4.4) 


Writing y in terms of x from the equation of the ellipse and 
integrating we get. 


V = 'n;b^(x ~ -^) 
3a'^ 


(4.5) 


The force per unit mass is therefore given by, 

se 
F 


■3: 

2 p 




5a -X 


where P is the density of the liquid. 


(4.6) 
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In a conservative system, force f is given by the expression, 

f = - (4.7) 

where 9^ is the scalar potential, 9^ is obtained by integrating 
equation (4.7) and is given by 


9- = - / P dx 


(4.8) 


Substituting P from equation (4.6) in equation (4.8) and 
integrating we get. 


, es^E ^2 

= - i 1 “ 1-¥-7 


( 4 . 9 ) 


4.2.2 Determination of pressure 


The effective pressure a P which is operative in 
restraining the fluid flow, is given by, 




(4. 10) 


where R^,R 2 are the radii of curvature of the ellipsoid and 
are given by the expressions (2.8) and (2,9) of Chapter 2. 

Y and r are the siirf ace tension and radius of metallic particle, 

Substituting y’ and y« (y' and y" = equation (4.10) 

and on simplification we get. 


w ( aA ) ( 1 /t ) ( r t ^ ^ 1 - 1 P )- r 3 

[(a/'b)2(1-U/a)2)+ ^ 


(4.11) 



^••2*3 Fluid flow equations 


For an incompressible liquid the continuity equation is 
given by, 

^ ’ q = 0 (4.12) 

where q is the velocity vector and is given by, 

q = i u + X V (4.13) 

where u and v are the components of velocity in x and y 
directions respectively. 

From the definition of circulation and Stokes theorem^ 

. * d.r = J^^n*(vx'q) (4.14) 

where indicates the integration over a closed contour c 
in the liquid and A the area enclosed by this contoxir. Since 
the liquid flow is irrotational, 

V X q =0 (4. 15) 

Hence from equation ,( 4. 14) we get, 

q • dr = 0 

98 

Equation (4.16) implies that. 


( 4 . 16 ) 
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where 9 is a single valued scalar velocity potential function. 
Substituting equation (4.17) in (4.12), the continuity equa- 
tion becomes. 


V^(p = 0 (4.18) 

For an inviscid liquid, the Havier-Stokes equation is given 

i= (4.19) 

where P is the body force per unit mass, the pressure 

gradient and P the density of liquid. 


In a conservative system, the body force is given by equation 
(4.7). Substituting equation (4.7) in (4.19)» 




( 4 . 20 ) 


For irrotational flow 
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the ( Dq/Dt ) can also be written as. 


Bi 

Dt 


V 


L2. . ± 
at ^ 2 


^q (4.21) 

The integration in equation (4.16) is evaluated at a fixed 
time. Hence, the time derivative can be written as. 


^c it ° (4.22) 

or 

^c ^Dt • + 9 • dq) = 0 (.4.25) 

Substituting ^ from equation (4.20) in equation (4.23) and 


integrating we get. 
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- 9^ ~ ^ + "i" ^ (4.24) 

Similarly substituting ^ from equation (4.21 ) in equation 
(4.23) and integrating we get, 


+ 4- q^ =0 


at 


+ 1 


(4.25) 


Equating (4.24) and (4.25) we get. 


3CP 

at 


) 


( 4 . 26 ) 


Substituting cp^ and P from equations (4.9) and (4.11) res- 
pectively and putting the dot product of q (equation 4.13) iu 

equation (4.26) we get, 

2 


at 4 P 




■Z 2 ^ 2 

3a -X 


i(u^+v^) - 




(4.27) 


In our analysis the closed contour is the boundary of the 
liquid droplet. From symmetry of the geometrical shape we 
consider only the first quadrant in the x-y plane, viz., 
o ^x £ a and o £ y £ b . 

As we have mentioned earlier, the boundary moves with time. 
Iherefore, one has to consider the dynamic boundary condition. 
This IS given by , 
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li i 
9t 1 , 

DO-undary 


[v-u- ® } 

boundary 


Initial conditions, at t = o, are 


u = V = 9 



(4.28) 


(4.29) 


According to our assumption, the initial shape is spherical. 

At t >0, from symmetry consideration, along the straight 
boundaries (x =0, y = o), = o and ^ = o respectively. 

Along the curved boundary , the functional value of 9 is obtained 
from the solution of the equation (4.27). Thus, the problem 
is to solve the laplace equation (4.18) consistent with the 
above boundary conditions. Since no analytical solution is 
available, we have solved the equation niunerically by using 
a finite difference scheme. 

4.3 Numerical Solution 

For ease of computation it is necessary to normalise 
the equations. The normalisation is carried out by setting. 


X = x/a 
Y = y/b 

and 

t* = t/t^ 

such that 

o < X < 1 , o < Y < 1 


(4.30) 

(4.31 ) 

(4-52) 

and o < t* <1 ( 4 . 33 } 
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By using expressions (4.30) and (4.31 ) the Laplace equation 
(4.18) becomes, 


, 2 2 

—J a_2. . 1_J. _ 0 

(a/b)^ 3X2 


(4.34) 


By using expressions (4.30) to (4.32) equation (4.27) becomes, 

,2 

In (— n) . |.(u2+y2) - 
3~X'^ 

^ [ (!■) > - 4 


I_ i£_ ^ 1 — o— 
7 3t* 4 » 


b^ 'b- 




■J 


where, u and v after normalisation become 


(4.35) 


u = 


V = 


1 

a 3X 


i 

b 3Y 


(4.36) 

(4.37) 


Using equations (4.30) and (4.31 ), the equation (4.28) becomes. 


b |~-+lT'Y=v-ur“ 


3t 


a 3X 


(4.38) 


Substituting expressions (4.30) and (4.31 ) in an equation of 
ellipse, we get, 


2 2 
X^ + Y^ = 1 


(4.39) 


Under this condition we must have, 

= o and evaluating from (4.39), the equation (4.38) 

becomes. 





X 
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(4.40) 

The above equation is valid for any value of o l_X 1 . There- 
fore, at X = o and using equation (4.32) equation (4.40) becomes 



^ =v| 

at* X=o 


(4.41 ) 


Since v is a function of time and its functional foim is not 
known it has been evaluated numerically. 

4 . 3.1 Kumerical scheme 

Figure 4.2 represents the first quadrant in the X-Y plane 
in which the equations have been solved numerically. The 
normalised X and Y aoces have been divided into 10 equal parts 
such that the steplengths are h = 0.1 and k = 0.1 respectively. 
Xo,Xo . . . , Xq ^ represent the points of intersection of 
the curve and the lines X = 0, X = 0,1, ..., X = 0.9. Similarly 
YqjYo . . . , Yq g represent the points of intersection of the 
curve and the lines Y =0, Y =0,1, ..., Y =0.9. These are 
the boundary points on the curve at which 9 's have to be 
evaluated at a panticular time so that the laplace equation 
can be solved at tha-t instant of time at the grid points as 
shown in figure 4.2. It is to be noted that the grid points 
marked with open circles near the boundary have unequal step- 
lengths. The following steps are involved in the numerical 
solution of the problem ^ . 
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Step 1 


Step 2 


Step 3 


: t^ is assigned an arbitrary value and the normalised 

time is divided intO' 10 parts such that step length 
p =0.1, t* =0 is set. 

t The first order partial derivative (of equation( 4.35) ). 
is repla.ced by backward difference. The finite 
difference formula becomes 


<p. * 


t +p 




+ t,p[| 


-r.2 

ee^E. * 
o t 


Iji ( - t(u^*+ v^*) - 


3-X 




(a^*/b^*)^( 1 )+( 1+X^) 


}- 



(4.42) 


^t*+p evaluated from above equation at the 
boundary points, ^q» •••» 9* 

E^* is given by the expression (2.18) of Chapter 2. 


: At every grid point (i,3) the Laplace equation 

(4.34) is replaced by the following finite difference 
formula, 




2? 


IjlI 


+ <p. 




(a^*/b^*)^ h^ 




(4.43) 


Two things to be noted in the above formula. 
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First, for any j, when i =0, we get 9 . 

1 f J 

9(-, . and 9 . in the first term of equation (4.43). 

From the houndary condition at t* > 0 on Y axis 
(|^ = 0), we get, 


^+ 1.0 - ^- 1.0 
23’^*h 


(4.44) 


From equation (4.44), the unknown 9 , . on the 
fictitious boiondary represented as dotted line in 
figure 4.2 is eliminated. Similarly, for any i when 

3=0,9. .is eliminated from the boxmdary condi- 

1 , “ 1 


tion at t* > 0, 


^ = 0 on X-axis. 


Second, for grid points represented by small 
circles in figure 4.2, the step lengths are xinequal 
and the situation around such a grid point is magni- 
fied in figure 4-3. The equation (4.43) is replaced 


2 ^ ^ 

7 ^'(X.+ 1) ^i-1,3 


Xi 


t 'Pi,3 + X. (x1h 

X ' 1' 1 


‘^i+ 1,3^ ^ 


^ X^+1 ‘^i,j-1 ^±,2 X^(x^+t ) 'Pi, 3+1 ^ ^ 


(4.45) 




Fig. 4.3 Magnified view of the surroundings 
of a grid point near the boundary for 
solution of fluid flow equations . 
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It is to be noted that for X. = X. =1 equations 
(4.45) and (4.43) become identical. Using equations 
(4.43) and (4.45) for all i,3's we get a set of 
algebraic equations which can be put in a matrix 
form. By using matrix iteration [Gauss iteration]^^ 
the 9 values at time (t*+-p) are evaluated, since 
the boimdary values evaluated in step 2 are valid 
for time t*+p. 

[in the present problem the matrix iteration 
scheme converges. This can be sho’tKi as follows. 

The absolute value of a diagonal element in any 
row of the matrix as stated earlier is given by, 


OK * 

I Diagonal element ! = — — 


.a.* p 
“t 


(f-) + 2k-2(n) 


( 4 . 46 ) 


The absolute sum of off --diagonal elements 


-2 T 

2h r 1 .r 1 Til 


II 




3 2 


2 h 


(f-) + 2k-^(-h) 




°t* 


^3 


(4.47) 
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step 4 


Near the hoiradary points either I or II or both 
terms will he transferred to the known vector side. 
Therefore, from equations (4.46) and (4.47) we get 
for the matrix. 


Absolute value of diagonal element 

^ Absolute sum of off-4iagonal elements (4-48) 

This (4-48) satisfies the condition for convergence 
of matrix iteration (Gauss-iteration)^^]. 

Velocity components at time t*+p are computed from 

<p. .'s as evaluated in Step 3 at time t*+p. At 

the interior points u. . and v. . are computed by 

i » 3’ 1 » 3 

the following finite difference formula, 


u. 




‘i+1,j “ ‘^i-1,3^ 

(4.49) 

'i,3+l " 

(4.50) 


At the interior points as represented by small open 

circles in figure 4.2, due to unequal step lengths, 

'■112 

the equations (4.49) and ( 4 . 50 ) are replaced by 
(see figure 4.3), 


u. 


t^ ^-(a^*h) f'* X (X +i) '?i-i , j 


i‘ i^ 




^i,3 ^ 


1 _ , ,1 

\TTxTTiT ^i+i ,3-* 


(4-51) 



step 5 


1 1 1 


»i, 3-1 - ^ ’■i.o ^ 


.+ 1 ) 


(4.52) 


It the boundary backward differences are employed to 

extrapolate the velocity components at points 

*■^"0 . 1 ’ • • • > g and Yq » Yq ^ ^ , . . . ^ Yq ^ g at time 

t*-t-p, from tp. . values evaluated in previous steps 
^ ♦ D 

(2 and 3). 

Replacing time derivative of equation (4.41 ) by 
backward difference we get, 


b . * 


= b .* + t.p V. 


t-_^p ^ 


(4.53) 


where v.*^_ I is the velocity component evaluated 
^ X=o 

at time t*+p at the point Xq in Step 4. As the volume 
of the ellipsoid remains constant, using as 

evaluated by equation (4.53), we get. 


3 2 

®'t*+p ~ ^ '^^t*+p 


(4.54) 

Since the potential drop between the two grains 
remains unchanged even after stretching, the field 
increases as a function of time. Prom figure 2,3 of 
Chapter 2 it can be shown that, at time t*+p, 

®t»+p = V V4+p> 

where s^*_^p = ■” 2(a^*_j_p ^ r) (4.56) 
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Step 6 
Step 7 
Step 8 
Step 9 

Step 10 
Step 11 


Step 12 


Step 13 


t* = t*+p is set. 

Steps 2 to 6 are repeated, till t = 0.9. 

The difference "between is^/2 - is evaluated. 

Depending on whether difference is positive or 
negative, time t^ is increased or decreased. 

t^ in Step 1 is replaced by t^ in Step 9. 

Steps 1 to 10 are repea.ted. This process is con- 
tinued xmtil we get two values of t^ such that 
in one case the difference (s^/S - is negative 

and in other case the difference is positive. 

Using the method of hisection t^ is evaluated and 
time t^ in Step 1 is replaced by this new value of 

7 - 

Steps 1 to 8 and then Step 12 are repeated till the 
difference as evaluated in Step 8 is within + 0.l2. 
and the time which is evaluated in Step 12 is the 
required time. 


The programme developed for the solution is given in 
Appendix 3. The programme is written in Fortran language with 
'some features permitted in DEC ^stem 1090 computer. 
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4-4 Experimental 

4.4.1 Preparation of glass and blowing of thin film 

A glass 'batch, of composition j 1 5 P 2 O 5 , 5 Bi 20 ^ (mole ^ ) 

has been melted. The melting is described in Section 3.5.1 
of Chapter 3- The blowing of thin films used for pulse response 
measurement has been described in Section 3.5.2 . 1 of Chapter 3. 

4 . 4.2 Pulse response measurement at room temperature 
4.4.2 . 1 Pulse circuit 

Figure 4.4 shows block diagram of circxiit configuration 
indicating the principle of pulse circuit operation, A standard 
pulse generator (Systronics Dual Channel, Model 1110-D) is used 
as a source of stimulus for providing a rectangular pulse of 
appropriate height and width to be impressed across the sample. 
However, the driving capacity in terms of power and current 
demanded by the sample of pulse generator is inadequate and 
therefore a buffer circuit combination consisting of Schmitt 
trigger interface and a solid state high speed (transistor) 
switch are interposed between the stimulus pulse generator and 
the sample. The regulated D.G. power supply (Aplab, Model 7152) 
is essential for providing the requisite supply voltages for 
proper biasing of Schmitt trigger and the power driver switch. 
In addition, it performs the essential function of determining 
the height or amplitude of final pulse that is applied across 
the sample. This is accomplished by varying the front panel DC 
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Block diagram of circuit configuration showing 
pulse circuit operation . 
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level control, provided on this equipment and monitoring it on 
the ORO which is also connected to the circuit at relevant 
mode for this purpose. The CRO made hy Hewlett Packard 
(Model 1725A| 275 MHz) is a dual trace oscilloscope with 
delayed time base and digital volt/ohm meter. 

Figure 4.5 shows the schematic representation of the 
circuit used for pulse response measurement. Transistor 
(ECH 055) is used as a shunt type of electronic switch and the 
current limiting resistor R^ ^ is in series with the parallel 
combination of transistor and the sample. The positive DC 
supply voltage (V ) is used for supply the pulsed current. 

When the input pulse to the base drive of turns it off 
from normally on -state, a pulse height of V. volts and 
duration of x seconds is impressed across the sample, conse- 
quently when the pulse height and width are correct the sample 
undergoes switching from the off to the on-state thereby 
conducting current of significant magnitude limited only by 
the series combination of resistor R^ ^ and the sample in the 
on-state. Finally, on termination of pulse dxiration x, the 
presently off -transistor turns on again but this time shares 
the current in parallel with the low impedance of the switched 
sample. To apply a pulse which first turns off on its 
leading edge and then turns it back to the on-state x seconds 
later, the requisite base drive current is supplied by the pre- 
driver stage consisting of transistor (2N 2904) which in 
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turn is driven by the composite combination of transistors 

(2N 708) and Q 2 (2F 708) and associated components connected 
in Schmitt trigger configuration. The Schmitt trigger circuit 
(Q| + Q 2 ) is a bistable regenerative swi'tciaing circuit designed 
to buffer the pulse generator from the load of pre -driver 
and driver stages (Q^ + Q^) and offers the additional important 
advantage of high speed switching action rith improvement in 
the time of transitions on the leading and. trailing edges. 

4 • 4 - 2 . 2 Measurement procedure 

For measurement we have placed the thin film sample 
vacuum deposited with gold on two faces between two metallic 
electrodes noused in a measuring cell as described in Section 
3. 3 . 6, 2 of Chapter 3. The rectangular pulse is provided by 
a standard dual pulse generator. The internal trigger is 
used to provide a direct output monitored on a high bandwidth 
delayed time base (275 MHz) dual trace CEO with sample dis- 
connected from the circuit. This is necessary to obtain a 
recurrent pulse train on the CEO screen which is held sta- 
tionary thereon enabling to adjust via pulse generator con- 
trols the width of the positive pulse position to the requisite 
magnitude prior to each switching operation. The precision and 
accuracy of pulse width measurement is enhanced by the availa- 
bility of the feature of time interval mode capability built 
into the specific CEO. The height of pulse act\xally driving 
the sample is Y , being determined by tlie DC regulated supply 




voltage, Next, mode of pulse generator is changed over to the 
manual trigger mode and the sample is connected to the circuit 
in shunt across the driver Q^. Now using the push hotton 
microswitch of pulse generator a single shot pulse is injected 
into the circuit causing the driver to apply a single shot 
pulse drive across the sample so as to subject it to switching 
action. Finally the sample is disconnected from the circuit 
and its resistance is measured on a low voltage ohm meter 
(1V DC excitation) so as to determine whether the sample has 
switched from the off -state to the low resistance on-state 
(memory state). First we fix the pulse height at 5 volts 
and sequentially increase the pulse duration from 10 micro- 
seconds to a few hundred microseconds by suitable interval. 

Each time after passing a single pulse the sample is checked 
whether it has switched or not. 

In the event of no switching taking place the pulse 
height is increased and the above sequence of experiments are 
repeated till the memory state is achieved. 

4.5 Results and Discussions 

Figure 4.6 represents the stretching of semi-major axis 

of the ellipsoid in molten condition by the critical switching 

field as a function of time for bimnuth metal of radius 

150 i and separation to diameter ratio 0.05. Fhe time required 

-1G 

to join up is 9.5x10 


second. It is seen from the figure 


Semi major axis 



g. 4.6 Stretching of two neighbouring elli 

as a function of time under the 

influence of an electric field. 
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"thati ini'tiall^ "the semi-inajor axis elongates slowly but 
later on at a rapid rate. This is because of the fact that as 
the distance between the two proximate grains decreases, the 
effective field increases (equation (4.55)) thereby causing a 
faster rate of deformation. 

Table 4.1 gives the switching time (tj^) as calculated 
from the solution of the time dependent heat equation in 
Chapter 3 and the time (t ) required for two proximate grains 
to coalesce as obtained from the solution of the fluid flow 
equations which has been described in this chapter for 
different diameters and separation to diameter ratios, 

¥e observe that t^ remains almost constant for a particular 
diameter even if the s^/d ratio increases. This arises due to 
the fact that a larger electric field is present in the latter 
case which lowered the time taken by the granules to stretch by 
the required amount. It is evident from the table 4.1 that t^ 
is less than tj^ by approximately four orders of magnitude. Thus 
the stretching dynamics is much faster than the kinetics of the 
heating of the sample due to Joule losses. This, therefore, 
justifies the assumption made in Chapter 3 for the computation 
of I-Y characteristics that once the required temperature pro- 
file is achieved within the sample, coalescence of proximate 
metallic granules takes place. 

Table 4.2 gives a few typical values of switching time 
obtained by pulse experiments. These are fommi to have values 
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Table 4.1 Calculated switching time for a glass of 
composition SOV^O^, 1 5P205>5Bi20^(mole^ ) 

for different d and (s^^/d) ratios 


Particle 

diameter 

i 

Separation to 
diameter ratio 
(s„/d) 

Switching time 
for sample of 
thickness 10 isn 

t^xlO sec. 

Switching time 
as calculated 
from stretching 
of grains 

Q 

tgXiO sec. 


0.02 

1.32 

0.95 


0.03 

0.72 

0.95 


0.04 

0.50 

0.94 

300 

0.05 

0.36 

0.94 


0,06 

0.28 

0.93 


0.07 

0.23 

0.93 


0.02 

7.39 

1.46 


0.03 

3.88 

1 .45 


0.04 

1.82 

1 .45 

400 

0.05 

1 .36 

1.44 


0.06 

1.07 

1 .43 


0.07 

0.89 

1 .43 


0.02 


2.03 


0.03 

5.01 

2.03 


0.04 

2.88 

2.02 

500 

0.05 

2.05 

2.01 


0 .06 

1.59 

2.01 


0.07 

1.29 

2.00 


0.02 


2.69 


0,03 

8.49 

2.68 


0.04 

4.15 

2.67 

600 

0.05 

3.41 

2.66 


0.06 

2.12 

2.65 


0.07 

1.72 

2.64 


700 


800 


900 


0.02 

- 

3.38 

0.03 

18,46 

3.36 

0.04 

5.86 

3.35 

0.05 

5.69 

3.33 

0.06 

2.72 

3.32 

0.07 

2.16 

3.31 

0.02 


4.13 

0,03 

- 

4.10 

0.04 

8.60 

4-10 

0.05 

4.82 

4.08 

0.06 

3.42 

4.06 

0.07 

2.65 

4.04 

0.02 


4.92 

0.03 


4.90 

0.04 

- 

4.88 

0.05 

6.36 

4.86 

0.06 

4.27 

4.85 

0.07 

3.23 

4.83 

0.02 

— 

5.76 

0.03 


5.74 

0.04 

- 

5.72 

0.05 

8.63 

5.69 

0.06 

5.31 

5-68 

0.07 

3.90 

5.65 


1000 
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around 50^-60 microseconds. The last two columns include the 
switching times t^^ and t^ respectively taken from columns 3 and 
4 of table 4.1. Columns 1 and 2 of table 4.2 give the switching 
fields as obtained by DC (from figure 3.6 of Chapter 3) and 
puls8 measurements respectively. In column 3 of table 4.2 

f 

the theoretically cs^lculated switching fields (taken from 
column 6 of table 2.1 of Chapter 2) are given. 

Comparing columns 4 and 5 of table 4.2 we see that the 
calculated t^ is in agreement within order of magnitude of the 
experimental switching time for particles of diameters 300 i 
and 400 i with separation to diameter ratios of 0.02 and 0.03. 

However » t is much smaller as compared to the experimental 
s 

value. Hence, switching at room temperature is controlled by 
the Joule heating step as we have mentioned earlier. However, 
it is possible in principle to get a switching speed of the 
order of nano -seconds provided the metallic particles are 

in the molten state. 

Comparing columns 1,2 and 3 of table 4.2 we see that 
experimental switching fields (DC and pulse) are within order 
of magnitude of the predicted fields for particles of diameters 
300 i and 400 i with separation to diameter ratios of 0.02, 

0.03 and 0.04. The differences in the switching fields as 
measured under DC and pulse conditions are attributed to the 
variation in micro structure depending on the melting and the 
blowing conditions of glasses. This agreement is reasonable 
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in view of the fact that a number of simplifying assumptions 
have been made in developing the model in Chapter 2. 

It is necessary to raise the temperature to the melting 
point of bismuth so tha,t the molten metallic bismuth granules 
can be stretched under the action of an electric field, it 
has been shown in Chapter 2 that the critical electric field 
for heating is less than the critical field for memory switching. 
Thus the electric field applied will raise the temperature 
level within the material to a value much higher than what 
is required for melting of bismuth granules. As the tempera- 
ture increases the glass structure becomes more open and also 
due to the fact that glasses are fairly low 

melting it is believed that the resistance offered by the 
matrix to the stretching of molten bismuth particles will be 
small as compared to the surface tension forces. We have, 
therefore, neglected any such opposing force by the gla.ss 
matrix to the granule deformation in our analysis. 



CHAPTER 5 


OPTICAI PROPERTIES OP 71NADITJI4 PHOSPHATE AID SILICATE 
GLASSES CONTAIIIHG BISMUTH GRAIULES 

5 . 1 Introduction 

Colloidal suspensions in glasses scatter light. Colloids 
of large particle size and sufficient concentration can 
generate turbidity in the glass which is undesirable. There- 
fore, to induce interesting optical properties highly absorbing 
materials are needed to predominate over the scattering. Metals 
being highly absorbing meet this condition provided the particle 
sizes are kept below 50 nm. Therefore, a proper choice of 
glass compositions and control of melting conditions can lead 
to the formation of ultrafine metal particles in a glassy 
matrix, 

Doremus^ ^ has studied the optical absorption of 
photosensitive alkali aluminoborosilicate glasses containing 
gold and silver. The latter are precipitated in the colloidal 
form by UV radiation followed by a suitable heat treatment. The 
spectral dependence of optical absorption has been analysed by 
lie theory. In the case of glass containing gold the corres- 
pondence between theory and experiment is good while in the 
other case the measured absorption compares reasonably well 
with that predicted by lie theory with regard to w>ne shape and 
position of the absorption peak. 
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Yellow, amber and red colours can be induced in the 
glasses by the process of staining. In this process compounds 
of silver or copper or a mixture of these are applied on the 
surface of the gla^ss. Then the coated glass is baked to a 
temperature in the range of annealing so that silver and copper 
undergo ion-exchange with the alkali metal-ions. The sub- 
sequent reduction treatment induces metallic colloids within 
the surface of the glass 'by nucleation and growth. The optical 
absorption of the borosilicate glass (composition similar to 

'Pyrex') stained with copper and a mixture of copper and 

1 1 5 

silver has been studied by Rawson Mie theory has been 

used to quantitatively account for the spectral trananission 
curve in the range of 0.4 to 0.8 ^im. The quantitative agree- 
ment between the calculated and the experimental transmission is 
good for glasses stained with copper alone, but the agreement 
is less satisfactory when a mixture of copper and silver is 
present. This deviation is attributed to the formation of 
copper-silver a3.1oy particles. 

Chakravorty et al. ^ ^ have studied the optical absorption 
of the float glass which has been subjected to ion-exchange with 
copper salt followed by reduction. The ion— exchange and 
reduction produce metallic copper within the surface of the 
glass. The characteristic absorption bands have been explained 
by Maxwell-Garnet t theory. Recently, optical studies in the 
visible range have been carried out on alkali containing 
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C -glasses and borosilicate glasses which have been subjected 

to ion-exchange in molten IglO^ followed by a reduction 

treatment in hydrogen gas at different temperatures and for 

.117 

various time periods . The absorption peak position has 
been explained on the basis of Maxwell— G-arnett effective 
mediim theory. 

The optical properties of cermet films produced by 

co-sputtering or evaporation (see Section 1.6 of Chapter 1) 

have also been studied extensively. The early measurement of 

optical properties of evaporated Au/SiO^ cermet films in the 

range 0.4 < X < o,7 has been reported by Hampe"*^®. Cohen 
119 ' 

et al have investigated in details the optical properties 
of CO -sputtered Ag/SiO^ and Au/Si02 cermet films in the entire 
range of 0.2 \m to 1.7 imi. lissberger et have 

studied the optical properties of MgT' 2 /^^ ajid Co/MgP 2 in the 
range of0.4< X < 0,6 tan. 

122 

Granqvist et al. have studied the optical properties 
of ultrafine gold particles produced by evaporation of pure 
gold in presence of air at a few torr. The percent trans- 
mission shows a characteristics transmittance dip at around 
0.6 lan. Extensive computer calculations have been carried out 
based on the effective medium theories of Maxwell-G-arnett as 
extended by Polder and van Santen (MC-PvS), of Himderi (HU) and 
of Bruggeman (BR) to explain the spectral dependence of optical 
absorption. It has been shown that in the limit of low 
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concentration of meta.! all the above effective medium theories 
predict the same result. By invoking the dipole-dipole inter- 
action between the metallic particles a better quantitative 
agreement has been achieved. The dipole-dipole interaction 
has been accounted for approximately by a set of effective 
depolarisation factors. Purther assuming that the micro structure 
consists of certain fractions of spheres, f.c.c. clusters and 
single strand chains quantitative agreement within a few percent 
has been obtained. In the case of ultrafine alximinium produced 
by evaporation in a gas mixture of helium and oxygen, the optical 
transmission characteristics have been explained within the 
MG formalism by suitably incorporating the dipole-dipole intera- 
ction and size dependent dielectric permittivity of metal par- 

!23 

tides coated with oxide. However the transmission dip occurring 
at 0.85 G is much less compared to the theoretically predicted 
one. This has tentatively been ascribed to defect scattering 
within the metal particles. Similarly, the optical properties 
of metal-dielectric^^’ composites formed by evaporation 
of metals - nickel, chromium, silver and gold have been inter- 
preted by using effective medium theories. 

Abeles and Gittleman^^^ have studied the insulator rich 
Ag/QiO^ cermet films and by applying MG theory they conclude 
that the latter is superior to BR theory to explain the data- 
However, Granqvist et al. have studied the applicability 
of MG~PvS, HU and BR theories on the experimental data of 
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119 

Cohen et al. . It has been foimd that no single effective 

medium theory can explain the experimental transmittance dip. 

However, this feature can atleast he explained qualitatively by 

incorporating the distribution of depolarisation factors in 

the MG type of effective medium theory. Optical properties of 

co-sputtered Cr^O^/Cr^^, EgO/Av?^ and co -evaporated 

have been investigated with a view to study the suitability of 

their application to photothermal conversion. Craighead 
34 

et al. have found a good correspondence between the prediction 
of MG theory and the optical absorption data of Ni/il^^ cermet 
films. Granqvist has carried out theoretical analysis based 
on MG and BR effective medium theories on the experimental 
optical data of MgO/Au^^’ Il20|^/Au^^^ and 
analysis shows that BR theory is superior to MG formalism in 
Au/MgO and Au/Al20^ cermet films but no such definite conclu- 
sion can be drawn in the case of MgP2/-^^ cermet systems. 
Granqvist et have carried out extensive computer 

calculations based on effective medium theories and applied to 
the case of silver in glassy matrix ’ and to the 

cermet system. 

Compared to an extensive investigation on the optical 
properties of metal -dielectric composite systems produced by 
sputtering or evaporation, such studies on glasses containing 
metallic granules are few. In this chapter we report on the 
characterisation of optical properties of oxide glasses 
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containing bismuth granules of dimensions of the order of a 
few hundred angstroms. Analyses of the type used in the case 
of cermet films have been applied extensively in the present 
investigation. 

5.2 Experimental Procedure 

5.2.1 Glass preparation and making of filmR 

The glass compositions which have been investigated for 
optical studies are given in table 5.1. 

The melting and preparations of thin films for first three 
glasses are similar to those discussed in Section 3.3.1 and 

3.3. 2.1 of chapter 3. Glasses 4,5 and 6 have been prepared 
from reagent grade chemicals. ®2®3 been added as H^BO^ 

Na 20 as Na 2 CC^, and SiOg and have been introduced as 

their oxides. The glass batch is accurately weighed, thoroughly 
mixed in acetone and melted in a hard fired alumina crucible 

in an electrically heated furnace in the temperature range of 
1300*^C to 1400^0. After homogenisation, the films of glasses 
are prepared by blowing as discussed in Section 3*3.2. 1 of 
Chapter 3. 

5.2.2 Reduction of thin films, 

Eigure 5.1 shows the reduction furnace schematically. A 
pyrex tube closed at both ends is inserted -through the 

electrically heated furnace. R represents 


mullite muffle of an 
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Tatle 5*1 Compositions of glasses 


Gla^s 

Fo- 



Composition in 

mole percent 


Vs 

^2^5 

SlOj 

B^O, 

2 3 

]?a20 

Bi20^ 

1 

80 

20 

- 

- 

- 


2 

80 

15 

- 

- 

- 

5 

3 

70 

20 

- 

- 

- 

10 

4 

- 

- 

64 

26 

10 

- 

5 

- 

- 

64 

18 

10 

8 

6 

~ 


64 

14 

10 

12 



Fig. 5.1 Reduction furnace . 
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a refractory "block on which the thin glass films S sandwitched 
"between perforated metallic sheets are placed. A thermocouple 
TC is placed near the sample and is connected to an Apla'b tem- 
perature controller. The reduction treatment has been carried 
out by passing hydrogen at a rate of 100 c.c./minute. Olasse 
no. 2 and 3 of the system have been reduced at 

200°C for 1*0 hour and 2.0 hours respectively whereas the glass 
no. 5 and 6 of the Na20-B20^-3i02-"Bi20^ system have been reduced 
at 300 °C for the duration of 1,0 hour and 1.5 hours respectively. 

5.2.3 Transmission electron microscopic studies 

The technique for electron microscopic sample preparation 
has been described in Section 3,3.5 of Chapter 3. The micro- 
structure and selected area diffraction (SAD) patterns of all 
virgin and reduced glasses containing bismuth have been taken 
by using Phillips EM 301 Electron Microscope. The SAD patterns 
have been analysed in a manner similar to that discussed in 
Section 3-3.5. 

5 . 2 . 3 . 1 Average Particle size determinati on 

Eor each sample a ntimber of electron micrographs have been 
taken. .A few hundred particles are measured. Then using the 
following formulae the average particle size and standard devia- 
tion have been calculated. 


X 


n 

S 

i=1 


X. n./F 
2 3 


(5.1) 
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I 

s = [ I (5.2) 

3=1 

where and are respectively the mean particle size and 
frequency of the ith interval, n and H are the number of 
intervals and the total number of particles respectively, s 
is the standard deviation. 

5. 2-3. 2 Calculation of volume fraction 

Assuming the precipitation of particles to be homogeneous 

throughout the glass matrix, we have used the point counting 

method to estimate the volume fraction (f^) by the following 
136 

expression , 



where n’ and U' are respectively the number of cross points 
coinciding with the particles and the total number of points 
present in the selected area of. the transparent graph paper. 
However, the volume fraction calculated by this method is 
greater than the actual one because of the superimposition 
of all the particles across the thickness on the two dimensional 
electron micro photograph. Assuming the particles to be 
spherical the actual volume fraction f of the particles can be 
related to f^ by the following expression , 

(5.4) 

em 

where t is the thickness of electron microscopic sample, 
em 
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5 . 2.4 Optical absorption stud i pa 

Poi* studying tiie optical absorption the tiiin glass filnis 
h.ave boon sand wit ciied bstween two iiai’d papers each having a 
square hole of about 5 mm x 5 mm punched in it. The optical 
spectra of virgin as well as reduced Y20^-”p205-Bi20^ glasses 
(glass no, 2 and 3 ) have been taken in the visible 2*ange using 
base glass no.1 as the reference, while for Ha20-B20^-Si02 -21203 
glasses (glass no, 5 and 6 ), base glass no. 4 has been used as 
reference. The spectra for the base glasses (glass no. 1 and 4 ) 
have been taken using air as the reference. All the spectra 
have been measured using Carry 17 D Spectrophotometer. 

5.2.4 . 1 ■ Calculation of optical absorption coefficient, a 

Carry 17 D plots the optical density (O.D, ) as a function 
of wavelength Prom this curve at any wavelength we get, 

(0,D.)gg^ple = ) observed ^reference ■ ^ 5 . 5 ) 

Prom Beer~lambert’s law which relates optical density to 
absorption coefficient a, we get, 

2 . 303 [( 0 »Pf )o]3gerved ^reference ^ (5.6) 

"'"sample 

where t represents the thickness of the sample, 

sample 

5.3 Results 


Plgures 5. 2 (a) and (b) show the plots of optical density 
as a function of wavelength in the visible range for vanadium 



Glass no. 
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phosphate base glass no.1 and silicate base glass no. 4 res- 
pectively . While glass no.l shows an absorption edge, the 
glass no . 4 is featureless. Figure 5*3 represents the plots of 
optical absorption coefficients as a function of wavelength for 
glass no, 2 and 3. Different reduction treatments applied to 
these glasses are indicated on the curves. It is observed that 
as the reduction time is increased for the same amount of bis- 
muth in the glass compositions the absorption coefficients 
increase . The absorption peaks for glass no«3 subjected to 
different heattreatments occur at the same wavelength 440 nm, 
while for glass no. 2 they occur in the range of 450-460 nm. 
Figure 5-4 represents, the plots of a values as a fimction of 
wavelength for glass no. 5 and 6 for different reduction treat- 
ments. The arrow head attached to the curve indicates the 
scale of a to the right. Like the vana,dium phospha.te system, 
it is observed that as the reduction time increases a values 
for both glass no .6 and glass no. 5 increase. Unlike vanadium 
phosphate system, each absorption curve shows two broad peaks 
in the range 500-530 nm and 420 nm to 430 nm respectively. 

Figures 5.5 to 5.16 show the electron micrographs of the 
vanadium phosphate and silicate systems under investigation 
for different reduction treatments. Seduction treatments are 
indicated along 'with the figures. Table 5.2 contains the micro 
structural details of "*"^0 silicate and vanadium phosphate 
glasses. It is observed from this table that volume fraction f 




Absorption coefficient Oi os a function of wavelength 
VoQr-.PtOc glosses containing bismuth. 


: Reduction temperature 



Fig. 5.4 Absorption coefficient 


Pig. 5.5 Electron micrograph, of a sample of 
glass no. 2 

X 117,000 



Pig. 5.6 Electron micrograph of a sample of 
glass no. 2 reduced at 200^0 for 
1 .0 hour 
X 84.000 







2.0 hours 
X 84 ,ooq__ 










glass no *5 reduced at 300*^0 for 1 hour 
X 117,000 ^ 




Pig. 5* 13 Electron micrograpb. of a sample of 
glass no. 5 reduced at 300*^0 for 
1 . 5 hours 
X 117,000 



jPig.5. l4 Electron micrograph of a sample of 
glass no. 6 

X 66,000 - 
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'512.5.15 Electron inicrograpii of a' sample of 
glass no. 6 reduced sl % 300^0 for 


1.0 hour 
X 117.00 0 



- T?-i CT R 16 Electron micrograph of a sample of 
.Pig. 5. 1b ^le^^ reduced at 300^0 for 

' 1.5 hours 

X 117,000 — - 


i 



Talale 5.2 Micro structural parameters for different 
reduction treatments 
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for gla,ss no. 2 and 5 varies in the range O.Ol to 0.04 and 
the average diameter varies in the range 75 to 155 angstroms. 
For the silicate glasses f varies in the range 0.02 to 0,05 
and the average particle size varies in the range of 80-160 
angstroms. 


Figures 5*17 and 5.18 show the SAD patterns of vanadium 
phosphate system for glass no, 2 and glass no. 5 reduced at 
200°C for 2 hours respectively. Figures 5.19 and 5.20 show the 
SAD patterns for glass no. 5 reduced at 300*^0 for 1.5 hours and 
glass no. 6 respectively. Tables 5.3 to 5.6 give the calculated 

interplanar spacings fi'om figures 5.17 to 5.20 res- 

1 0R 

pectively. Standard values of metallic bismuth have 

been included in each table for comparison. In 

system (figures 5.5 ■- 5.10) globular structure is more prominent 

in figures 5.8 and 5.9, while in silicate system (figures 

5.11 - 5.16) globular structure is more pronounced except in 

figures 5.15 and 5.16. However, in all figures certain fractions 

of interconnected chains are discernible. 


5.4 Discussions 

Before discussing experimental results as given in the 
last section we briefly describe the salient features of the 
different theories available for explaining the optical behaviour 
of inhomogeneous materials. 



Fig. 5. 17 Electron diffraction pattern of 
sample of glass no, 2 




, Pig. 5. 18 Electron diffraction pattern of a 
sample of glass no .3 reduced at 
200OC for 2 hours 
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Fig. 5- 19 Electron diffraction pattern of a 
sample of glass no. 5 reduced at 
300^0 for 1 . 5 liours 



Fig. 5. 20 Electron diffraction pattern of a 
sample of glass no. 6 
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Table 5*3 


Comparison of calculated d, , , 

hkl 

glass no. 2 with standard d^^^^ 


■values of 
value s 


Calculated d 

i 


hkl 


Standard d, , , for bismuth 

s. ^ 


108 


2.20 

2.27 

1 .86 

00 

1.31 

1.33 

1-11 

1.12 

0.83 



Table 5.4 Comparison of calculated dj^jj.3_ values of 

glass no. 3 reduced at 200°C for 2.0 hours 
with standard values 


Calculated d 

i 


hkl 


Standard d ^^ -^ for bisnuth 

i 


108 


2.17 

2,27 

1 .84 

1 .87 

1.30 

1.33 

1-10 

1.10 

0.81 

- 
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Table 5.5 Comparison of calculated values of 

glass no. 5 reduced at 300°C for 1,5 
with, standard values 


Calculated 

i 


,.108 

Standard for bismuth 

. 1 

2-20 


2.27 

1 .86 


1.87 

1.53 


1.33 

1 .12 


1.12 

Table 5.6 

Comparison 
glass no. 6 

of calculated d^^-^ values of 
with standard 

Calcula.ted 

i 


. . 4 - 1,108 

Standard for bismuth 

i 


2.20 
1 .89 
1..33 


2.27 

1-87 

1.53 


1.12 


1.12 
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5-4*1 Mie Theory 

137 

Mie equations have been derived by applying Maxwell’s 
electromagnetic equations to small conducting particles whi<±L 
absorb and scatter light, the dimensions of the particles being 
smaller than the wavelength of the optical radiation. For 
identical spherical particles the solution of electromagnetic 
equations for each particle is given by 

2 

Qgcat = 5 ^5.7) 

m +2 

2 2 

Q , = Imaginary part of 4A(%^) + fr * 

m^+2 m'^+2 


2it:rN 

where A = — r- ^ 


^ m'^4-27m^+38 j 

2m^+3 


(5.8) 

(5.9) 


r being the radius of particles, I the refractive index of 

O 

glass and the wavelength. 


m 


I-iK 



(5-10) 


with K and K being the real and imaginary parts of optical 
constants of the particles. and are respectively 

defined as the proportions of light being scattered and absorbed 
by each particle divided by the cross-sectional area (-nr ). 

The optical density can be related by the following expression. 
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OD 


1 

4 


(5.11) 


where “ *^scat '^ahs* ^ the weight fraction of the 

metal particles per imit volume of the composite, p is the 
density of the metallic particles and 1 is the sample thickness. 
Since QgQg_-{; is related to the 4th power of radius of particle 
and the primary term of is linearly dependent on the parti- 

cle radius for sufficiently small particle radii (aroimd 20 nm) 
the contribution of will be dominant. In this small 

particle size range, as can be seen from equations (5.9) and 
(5.11) the optical density is independent of particle size. On 
the other hand the contribution from becomes significant 

when the particle size increases. Hence for large values of 
particle radii the OD is dependent on the particle size. In 
the latter case it is likely that the particles will deviate 
from spherical shape. Mie theory has been modified by 
Gans^^^’^'^^ to take into accomt the effect of deviation from 
sphericity, 

5.4.2 Effective medium theories 


Effective medium theories have been used quite extensively 
to metal-dielectric microparticulate systems as we have dis- 
cussed in Section 5.1. 

The simplest effective medixmi theory is the one which has 

lit 

been developed by Maxwell-Gamett‘ . The effective permeability 
has been derived on the assumptions that particles are identica.! 
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spheres, they are sufficiently separated from each other to 
ensure independent scattering and the Lorentz correction is 

applicable. Under these assumptions the effective permeability 

MGx . . , 

(e ) IS given by, 


e 


ja. 


£ - e 


+ 2s 


m 


e +2£ 


m 


m 


(5.12) 


where is the dielectric permittivity of the matrix, f the 
fill factor and e the permittivity of an individual particle. 


1 42 

It has been shown by Onsager that only part of the 

field is effective in directing the dipoles in a dielectric 

1 45 

matrix. Polder and van Santen have incorporated Onsager ’s 

t f 

reaction field to derive an improved effective medium permea- 
bility for particles of ellipsoidal shapes. Combining their 
calculations with MG approach the following effective medium 
permeability can be obtained^ 


^IG-PvS 
£ = £ 


IE 


1 + 2/3 
1 - 1/3 r 


( 5 . 13 ) 


where f . is the fill factor of the particles belonging to the 

jth class, -such that X-f- = f, where f is the fill factor. 

6 . is related to the polarisability of the particles in jth 
3 

interval and for ellipsoidal shape is given by 

3 


122 


, e . - e 

6 , = i I 


k=1 




( 5 . 14 ) 
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where is the size dependent dielectric permittiirity of 
the particle, e is the effective permeability which can be 
calculated from any effective medium theory. Lj^’s are the 
triplet of depolarisation factors. It can be seen that for 
equal-sized spheres with = l/3, the equation 

(5-l3) reduces to the Maxwell-Sarnett expression (5-12). Both 
MG and MG-PvS theories are applicable for small values of the 
fill factor. 


Another effective medium theory has been put forward by 
1 44 

Bruggeman . In this a particle is considered to be embedded 
in an effective medium and its properties are determined self- 
consistent ly . To achieve this one has to solve for local field 
around the particles such that the fluctuations average out to 
zero. Prom these requirements the effective medium permeability 
is given by, 


1-f + J- s.f.6. 

JBR ^ ^ 5 Ll-L 

m o 

1 - f - r -6 . 

' 3 3 0 3 


(5.15) 


where 6. is given by the expression (5.14). In principle this 
3 

theory is applicable to any fill factor, as it gives equal 
weightage to the components of the inhomogeneous medium. 


Hunderi^'^^ has derived an expression for the effective 
peimeability from the back scattered field and it is given 
by, 
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" = ^ 




’ - 4 '//a 


(5.16) 


6^ being given by ihe expression (5»14)« Tiie iiieory is 
applicable for small fill factor. 


Recently both MG- and BR theories have been derived by 
146 

Riklasson et al. from an unified approach by specifying two 
types of random unit cells which represent two types of micro- 
structures and on the requirement that these cells when embedded 
in an effective medium will not be detected by an electromagnetic 
wave . 


In the small particle size range, the grain boundary 

scattering should influence the mean free path of conduction 

electrons and hence the dielectric permeability will be 

particle size dependent. Secondly when the particles are not 

far apart the dipole-dipole interaction cannot be neglected. 

This interaction can be included by replacing the 6.'s in 

equation (5.14) by a set of equivalent depolarisation factors 

defined for different geometrical configurations of identical 
1 22 

spheres 

5.4.3 Validity of effective medium formalism 

The criterion for the validity of effective medium model 

1 P9 

as examined by Granqvist' is that the results predicted 
should he independent of experimental condition used. When a 
plane wave with aa amplitude is incident on a thick slab 
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containing independently scattering particles the hack 
scattered field has been given by Van de Hulst^'^'^ as 
iTcUoPpSCTt )/k where is the particle niunber density^ S(tc) 
the scattering amplitude function at an angle % and k is the 
wave vector. The refractive index of the medium corresponding 
to this field is given by^'^^, 

% 1 1 - i2x PpS(Ti)/k^ (5.17) 

In the transmission mea.surement the refra-ctive index is giTen 

hy y 

1. 1 ■> 12 % PpS(o)A^ (5.18) 

The ratio of S(n:)/S(o) as calculated by 7an de Hulst^^"^ ■ 
from the Mie theory is given by, 

S(ti;)/S(o) = 1 + o(k^x^) (5.19) 

where e is the dielectric permittivity of the particles. In 
order for the effective medium approach to be valid, therefore, 

and should not differ appreciably from each other. This 

K i 

means that S(ii;)/S(o) should be around 1. Tor the real and 

imaginary parts of the second term of equation (5.19) to be 

smaller than around 5 percent, the particle size will be of the 

1 22 

order of 10 nm or less for gold . The 3rd row of table 5.2 
shows that the particle size is of the same order of magnitude 
in the present investigation. Therefore in the following 
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sections we have used effective medium theories to interpret 
the optical absorption data as reported in Section 5.3- 

5 • 5 Interpretation of Experimental Results 

From tables 5.3 to 5.6 we see that the calculated 

values compare well with the standard d, , , values of metallic 

JlKl 

bismuth. The SAD data thus confirm that the particles observed 
in the different glass matrices consist of metallic bismuth. 

5*5.1 Computation procedure 

The frequency dependent dielectric permittivity of bianuth 
is of prime importance in interpreting the optical absorption 
data by effective medium theories. To start with the frequenc3r 
dependent bulk dielectric permittivity of bismuth has been 
calculated from the index of refraction and extinction coeffi- 
cient given as a function of wavelength in reference ([148]). 
The size dependent dielectric permittivity of blanuth has been 
obtained from the bulk data by incorporating the size dependent 

electron scattering in the Drude free electron part of the 

1 22 

dielectric permittivity as follows 





Drude 

^exp 



+ 


Drude/ \ 


(5.19a) 


where « is the angular frequency of radiation used. The last 
two terms in the right hand side of the above equation are given 


^ 7 , 
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Drude/ \ . 2 , , , , 

Sxp (“) = 1 - + lAi,) (5.20) 

^Drudej^j _ ^ (5-21) 

where and are the bulk plasma frequency and the mean 

electron life time respectively, m is the apparent plasma 

px. 

frequency of particle of size x. and -t . is given by, 

3 3 


tJ’ = + 2Tjyx. (5.22) 

where Vj,.j^ is the bulk .Fermi velocity of electrons and x^ the 

particle size of the jth class. A value of = l6.0 e¥ 

has been taken from reference ([149]). Values of = 2.838eV 

8 

and Vp.|^ = 1 .822x10 cm/sec. have been calculated using free 
1 50 

electron model . 


Separating the real and imaginary parts of the size 

dependent dielectric permittivity (e. = e..’ + ie^**) from 

J J a 

equations (5.l9a-)to (5.2l) we get after simplification, 


e . ' ( w ) 
3 ^ ^ 


e' (m) + i 
exp ^ ..2 


2 


(I) 


MX 


+ iAi3 “ 3 


(5.23) 


Sj" (“) 


05^ .[(x. + 2i:, Vt^)/t. X . 3 
u / \ . bxi‘-'‘ n b Fb^' b .t* 

b „ o . . . 2- 


exp 


a,[a)A{(x.+2T^Vj.^)A^X^ rj 


‘^^+{iA|) 


(5.24) 
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5. 5*1*1 MG effective medium theory 

Simplifying equation (5*12) and then separating the 

MG 

•1 + ^^2 


. . , „ — MG / -MG -MG -MG v j. 

real and imaginary parts of e (r = + ie“ ) we getj 


-MG 

^1 - n2 


AC + BD 


■where , 


C^- + D‘ 


-MG _ BC - AB 

Ep - 2 2 


A = e^(l+2f)£' + 2£^(l-f) 


B = ( 1 +2 f ) £ * 


C = (l-f)£' + (2+f) 


D 


(l-f)e'’ 


(5*25) 

( 5 * 26 ) 

(5.27a) 

(5.27h) 

(5.27c) 

(5.27d) 


In equations (5. 27), £' and £» are given by equations (5*23) 
and (5.24) respectively in which is replaced by the average 
particle size as computed by equation (5*1)* The frequency 
dependent and have been computed using equations (5-25) 

to (5.27). Bor bismuth in silicate and in ^ 20^-5205 systems 

4. 1 „ o o[l 28 ] 

the dielectric permittivity of the media are taken as 2. 
and 10 ^^ respectively. 

5 . 5 . 1 . 2 MG-PvS and BR theories 

Equations (5.15) and (5.15) contain size dependent fill 
factor which can he related to f hy the following equation. 



(5.28) 


^3 = "d" 


where 


^3 weight factor and is given by, 

3 

X. n . 

W . = — ■ ^" ■■ 

^ x^n. 

3 3 3 


(5-29) 


where n. is the n-umher of particles in jth class. 

From equation (5.13) separating the real and imaginary parts 
of the complex quantities, we get the real and imaginary parts 
of the effective medium permeability in the M6-FvS theory as 
f o Hows ; 


where , 


6 '. and 

J 




= 8 


+ A^Dij 

“ 0^ + 


A, = 1 + 5 f 


B, = ff i.w.s:- 


'1 


, - i f r , 
J ^ 5 " 3^3 


6* of the above equations 
3 

5 EGi, + 

aL 




(5.30) 

(5.31) 

(5.32a) 

(5.32b) 

(5.32c) 

(5.32d) 

(5.32) are given by. 


(5.33) 
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where , 








(5.34) 

(5.35a} 



(5.35b) 


+ 




(5.35c) 


% = Vj + (5.35d) 

In the above equations e’. and are given by expressions 
(5.23) and (5.24) respectively. In the Briiggeman formalism, 
the expressions (5.30) to (5.35) will be identical with 
■4^ ,0*1 and of equations (5.32(a) to 5.32(d)) being properly 

replaced by using the Bruggeman effective medium equation 
(5.15). 

In order to compute the value of and Eg from equations 
(5.30) and (5.31), we notice that the equations are of the trans- 
cendental type as they involve e^ and ^2 rrght hand side 

through 6'. and 6'.' respectively. So e". and ^2 have been found 
3 3 

iteratively as follows. Pirst in the right hand side is 
assigned the value of and we set £2=0. Then e| and £2 
are calculated from the equations (5*30) and (5.31 )* and £2 
of the right hand side of equations (5.30) and (5.31) are 
replaced by the new values as calculated in the previous 
iteration and e^ and £2 are calculated from equations (5.30) 
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and (5«3l)-. This process is continued until the difference 
between two successive iterations is less than a preassigned 
error. In our calculation we have used an error of 0.1. l|^’s, 

the depolarisation factors for sphere are taken from the 
reference ([122]). Once the frequency dependent and ^2 are 
computed, the optical absorption coefficient (a) is then given 

a(a)) = (u)/c)(e2/E^^) (5.36) 

The frequency dependent a values as computed as a function of 
wavelength X from MG, MG-PvS and BR theories for the '^2^5”^2^5 
and silicate systems under investigation are plotted in figures 
5.21 to 5.32. In each figure is included the corresponding 
experimentally determined a values for comparison. Table 5.7 
contains the theoretical and experimental a values corresponding 
to peak positions. It is observed that computed a values by 
BR and MG-PvS theories are almost identical as are evident from 
figures 5.21 to 5.32. So in table 5.7 we have only included a 
values computed by MG and BR theories. The theoretically com- 
puted a values for vanadium phosphate glasses containing bis™ 
muth show a small absorption peak at 650 nm while in the case 
of silicate glasses containing bismuth the calculated ot values 
show a small dip. However, no peak is observed experimentally 
in either of these systems around this wavelength. 
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Fig. 5 23 Comparison of experimental and thcorctico 
absorption coefficient a for gloss no. 2 . 





Fig. 5.24 Comparison of experimental and thcoretica 
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Fig. 5.32 Comparison of experimental and theoretical absorption 
coefficient (X for glass no. 6. 


Table 5.7 Comparison of experimental and theoretical a values for different 
glasses 
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In system as the reduction time increases the 

experimental a values at peak positions increase {2nd row of 
table 5»7 and figure 5.3). This is because of the fact that 
with, increasing reduction time the amount of precipitated 
bismuth becomes larger. It is evident from the second row of 
table 5*2 that the volume fraction (f) increases wilii an 
increase in reduction time. From the 5th and 7th rows of 
table 5.7 we see that a values as calculated by MG theory 
show two peak positions at 450 nm and 410 nm while BE theory 
predicts one peak at 460 nm. Beyond 440 nm the values of 
as calculated by BE and MG~PvS theories become negative. 
Experimentally we could not observe the second peak as the 
experimental OD could not be measured because of the full 
opening of the slit in Carry 172 beyond about 410 nm. Com- 
paring the 1st and 5th rows of table 5.7 we see that the 
first peak position as predicted by MG theory agrees with the 
experimentally observed value within + 10 nm. The peak 
position predicted by BE theory (9th row of table 5.7) agrees 
with the experimental one within an error of 20 nm. MG theory 
predicts results which are in better agreement with experiment 
than those obtained from BE and MG-EvS theories. This is 
evident from figures 5.21 to 5.26. It is observed that calcu- 
lated a values by MG, BE and MG-PvS theories increase with 
increase in reduction time except for glass no. 2 reduced at 
200°C for 2.0 hours. From 2nd row of table 5.2 we see that f 
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increases with increase in reduction time except for the one 
just mentioned. Except for this sample the theoretically 
observed trend of a values with reduction time is in agreement 
with experimental results. 

In the silicate system we find two broad absorption 
peaks* This is evident from figure 5*4* Prom 1st and 3rd ^ 
rows of tables 5*7 we see that two absorption peaks occ\nr 
in the range of 500—530 nm and 420 to 430 nm respectively. 

The corresponding calculated peak positions by MG theory occur 
at 470 nm and 430 nm respectively (5th and 7th rows of table 
5*7) and those computed by BB. theory occur at 530 nm and 460 nm 
respectively (9th and 11th rows of table 5*7). Therefore, it 
is clear that no effective medium theory by itself can predict 
the peak position accurately. 

Prom the study of table 5.7 and also from figures 5-21 to 
5.32 we find that the discrepancy between Ihe theory and the 
experimental data is more in silicate system than in the 
7^0^-P„0r- system. The calculated a values in the silicate 
system (6th, 8th, lOth and 12th rows of table 5.7) are an 
order of magnitude greater than those observed experimentally 
near the peak positions (2nd and 4th rows of table 5.7). This 
is believed to arise from the possibility at an over-estimation 
of the fill factor f as determined from the electron micrographs 
of the specimens concerned. The dark phase in the micrographs 
may consist of a metal phase and biaauth oxide rich glass 
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phased 52^ However, an accurate delineation of such micro- 
structural features is not possible by the experimental 
technique employed in the present investigation. 

Experimentally we observe from the 2nd and 4th rows of 
table 5.7 that a values increase with the increase in reduction 
time in silicate system. However, for glass no. 6 it is found 
that the calculated values show an opposite trend (6th, 8th, 
lOth and 12th rows of table 5.7). Ihis anomalous behaviour in 
predicted a values, though not magnitudinally very different, 
for the glass no. 6 and glass no. 2 as mentioned earlier can be 
explained as follows. Figure 5.33 gives the plots of a values 
as calculated by MG theory for different particle sizes in the 
range of 80 to l6oi with a fixed fill factor of value 0.02, 
Figure 5.34 shows the plots of a values as calculated by MG 
theory for a fixed particle size with different fill factors 
in the vanadium phosphate as well as in silicate systems. We 
see from the these two figures that fill factor has a more 
pronounced effect on the calculated a values than the particle 
size. The fill factor f “ (equation 5.4) where 

t is the thickness of the specimen under observation in the 

em 

transmission electron microscope. A value of lOOoX has been 

assxiDioci for ij wiiich. is, uppor ^ limif of elocfroB fraxis-^ 
oin 

mission. However, it is quite likely that it may be less than 
lOOOi. In that case the calculated f will be less than the 
true value and as a result there will be an underestimation of 


the a values. 



ou ssDiq 



etical plots of optical absorption coefficients OL 
ength X for different particle diameters (calcul 
G Theory) . 
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Theoretical plots of optical absorption coefficient a 
vs v^avelength X for efifferent values of fill factors 
(calculated by MG Theory! 
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So far we have assumed that the metallic phase consists 

either of spherical granules of some average particle size 

or of a distribution of spherical particles (BE and MG-PvS 

theories)* However, the inspection of micro structures 

(figures 5*5 to 5.16) suggests that elongated interconnected 

particles are also present. Therefore, the assumption that 

the metallic . bismuth phase consists of g fraction of 

spheres, c fraction of single strand chains and x( x=1”5-?) 

fraction of f.c.c. clusters is reasonable. This consideration 

effectively takes into account the dipole-dipole interactions 

between the particles as they are closely situated. This is 

accomplished by replacing the I'jj.'s in equation (5-14) by 

effective depolarisation factors corresponding to the above 

122 

mentioned geometrical configurations . It has been shown by 

1 22 

Granqvist et al. that there exists an optimum combination 
of g , g to obtain the best agreement between experimental and 
. calculated a values subject to the following constraints, 


5 >0 

(5-37) 


(5.38) 

c < t 

(5.39) 


Thus the problem takes the form of a constrained optimisation. 

, I ■ 

¥e have chosen the MG-PvS formalism and glass no. 2 reduced at 
200^0 for 2,0 hours for such an analysis. 
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5*5*2 of co n strained optimisation and nmaerical 

computation 

The objective function to be minimised subjected to the 
constraints specified in equations (5.57 ~ 5,39) is given by, 

n 

f(l) = I _ ^expj2 


m 


(A)=1 


©XT) 

where and are the theoretical a values as calcu- 

lated by MG-— PvS theory and the experimentally observed a values 
respectively at an angular frequency a and n being the number 
of frequencies at which the a value is calculated. X is a 
vector consisting of C - the optimisation variables. The 
above constrained optimisation is converted into an equivalent 

•4" 

unconstrained optimisation by augmenting the function f(X) 
by a penalty term within the frame work of internal penalty 
function method The equivalent function 9(J, r^) is 

given by. 


m 


(p(X,rj^) = f(i) - r^ I 


j=1 




(5.41 ) 


where r^^ is the penalty parameter^ gj(X) 5 . 0 is the constraint. 
In our case the last teim of equation (5.41 ) is given by, 


m 

I 

j=1 


± 




_ r » i « i 

- L r t 


+ 


t 


S + C-1 




(5.42) 



186 


For evaluating the hy MG—PvS theory and will 

be given by equations (5.30) and (5.3l) with and 0^ 

of equations (5.32) replaced by the following expressions to 
take into account the ^ fraction of sphere, c fraction of 
single strand chain (s.c.) and X fraction of f.c.c. 
cluster (f.c,). 


A. = 1 + I- f[S I w.6^ I +^6't + ] (5.43a) 

‘ ^ J 3 jgpiiepe |s.c, ff.c,-' 

3 


I W. 6! 


+ + X6‘M ] (5.43b) 

sphere Is.c. if.c. 


0, = 1 - ^ f[5 I 


+ C6M + X6' 
sphere Is.-C. 


P^c.K5.43c) 




(5.43d) 


where 


^3 


and 6' I a^® respectively the real 

s.c. jf.c. 


( sphere . - 

parts of tue polarlsabilitjr corresponding to spheres.olagle 

strand chains and f.c.c. ctasters. Similarly, 

6” I and 5" U . represent the imaginary parts of the 

pohl-sahility ;c;responding to the geometrical oonfi^ationa 

• strand and toe f.#!.P.clusters. 6. and 6^ 

spheres, single strand cnauo 

'nnc, f5 33) to (5.35) by introducing \ b 
are given by expressions (5.3:^; 

,122 The 

for sphere • 
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corresponding to different geometrical configirrations have 
been computed, as we have mentioned in Section 5.4.2, t ha t 
these are made up of identical spheres’ Therefore, in the 
case of single strand chains and f.c.c. clusters an average 
particle size has been used to compute the values of 6' and 6". 
After evaluating and n-^erically, the a values are computed 
from equation (5.36). ^ 


5 • 5 . 2 . 1 Numerical scheme 


The following iterative scheme is used to optimise the 
function 9 ( 1 , r^^). 


Step 1 


Step 2 


Step 3 


Step 4 


The computation is started with three initial 
feasible values of X, ,]£„ and X, such that the 
strict inequality constraint is satisfied. In 
initial value of r^ > o is taken, k = 1 is set. 
9 (X,rj^) is minimised using Simplex optimisation 
technique’ Whether the minimised X^* obtained 
corresponds to the minimum of the original problem 
or not (f(X)) is tested. If so the process is ter- 
minated. Otherwise the next step is followed. 

New penalty parameter is computed by 

T = cr, , where c < 1 . 

^k+r h" ^ 

With the minimum value of X^^* and setting k = k+1 
step 2 is repeated. 
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The computer programme in fortran language is given in 
Appendix 4 . 

Fig^ure 5.35 shows the plot of a values as a function of 
wavelength ^ for an optimum comhination of 26 percent spheres, 
0.6 percent single strand chain and 73.4 percent f.c.c. clusters 
within an error of + 4.0 percent for glass no , 2 reduced at 
200°C for 2 hours,- In this figure sere also included the 
computed values from MG, and MG-PvS theories and the experimen- 
tal results. Figure 5.35 shows that a hetter agreement is 
obtained between theory and experiment when the micro structure 
is assumed to consist of a combination of spheres, single strand 
chains and f.c.c. clusters rather than a simple distribution of 
spheres- It may be mentioned that the above optimisation 
procedxxre when applied to the silicate glass system does not 
bring about any significant improvoaent in the calculated a 
values vis-a-vis the experimental data. This is attributed to 
the possibility of an overestimation of the fill factor as 


discussed earlier. 
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Fig. 5.35 Comparison of experimental absorption co 
icient a with theoretical values computed 
by different models for glass no. 2 reduc 
at 200® C for 2:0 hours. 




CHAPTER 6 


COFCIUSIONS 


Based on electrical and optical studies of vanadium 
phosphate as well as silicate glasses containing bismuth the 
following conclusions are drawn. 

( 1 ) A particle stretching model in which molten metallic 
granules stretch under the influence of an electric field 
can explain semi-quantitatively most of the characteristic 
features of memory switching in oxide glasses containing 
metallic bismuth granules. 

(2) Depending on the sample thickness the vanadium phosphate 
glasses containing bismuth can ahow either memory or threshold 
switching behaviour. 

(3) The computed current '■voltage characteristics especially 
in the negative resistance region compare reasonably well with 
the experimental results for vanadium phosphate glasses 
containing bismuth. 

( 4 ) The switching time as computed from the numerical solution 
of time dependent heat equation in the presence of Joule losses 
agrees well with experimental switching time obtained by pulse 


measurements. 
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(5) switching time as computed from the numerical, 
solution of time dependent stretching of molten metallic bis- 
muth particles in vanadiimn phosphate system, is found to be 

very small as compared to the experimentally determined switchisag 
time from pulse measurements. However, the fast kinetics of 
stretching indicates that if the metallic bisnuth particles are 
in a molten state switching time of the order of nanoseconds 
can be achieved. 

(6) MG-, ,MG--PvS and BE theories predict the optical absorption 
peak positions which are in fair agreement with the experimentally 
observed values for both silicate and vanaditm phosphate glass 
systems containing bismuth granules. 

(7) mg- theory predicts results with regard to the absorption 
coefficient a which are in better agreement with experiment than 
those obtained from MG-PvS or BE theory in vanadium phosphate 
glass-*ffletal system. 

(8) A typical computation using constrained optimisation 
technique indicates that better quantitative agreement (within 
+ 4 percent) with experiment over the entire spectral range can 
be achieved by assuming that the vanadium phosphate glass-metal 
composites consist of certain fractions of spheres, f.c.c. 

clusters and single strand chains. The microstructural 
features of most samples tend to corroborate this model. 
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r!Mt£RPOLATTUN*nF FTHST"tJPDER DERIVATIVE OF TEMPERATURE 
BY RFGULA FAGsr MfcTHOn 


TPTs:(TPp*UlFM-TPNtDIFP)/(OIFN-DIFP) 

r:<\r,G gopes 

TF([)IFl..l,E,y,l)GO TU 74 
G'i tn 16 . ^ 

TFCOTFT.GT.O.) TP1=TPI*.9 
RETURN 
K .-i D 

RUPHOUTIME GOPES 


niMENGTUM E(2b) ,TPC200) ,T(2001 
r,j'4Mnr*/At/F,TriIC,K 
rii‘^Mnrj/A2/TTfFUN „ 

ri}^^MnN/A3/TPI,TMB,Tl.TMX»N0 

CU4WOfi/A4/DlF,DlFl TPP,TP2,Dp;P-01EW 
riJMMON/AS/DlFT^DTFTl ,THlCPf THlCN^Dlf'TP.OIl TU 
S^NDsMD 
0=1, /HMD 
0 ,^i=Mn+i 
Tprn=Tpr 
T C 1 ) = T I 
DO 60 I.=7. »MN 

TJ(TTi‘GT!50,3TPl*TPl + l,E2 
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60 


'?U 

?b 


40 
45 
1 35 
700 


t05 
7 3 


CAM, H-ISK 
OHl=f„/^,)4!{.Tf),( , 

OKr. = nMrP(6-l)+BBl) 

T’lMTru-l. I +DKT, 

^• CO TO 30 

mAjiIj Kn^r. / ■ 

n.)=tH/2. )*PUN 

AK.S=!./3.^(AAt4.<=BBl + Dn) 

T,U’? T * i'iT^ ^ (TP(T,-1 ) +AK ) 

rtMb)=TP(L-nf AKS 
CU'U'TOM*: 

OiF=TtMN)-T( I ) 

’UF1=ARS(0IF) 

TK(UIF,GT,0,)GO TO 20 
Tp'l=tPC 1 ) 

OtFr4 = i>TF 

Cij rn 25 
T,>t>sri>(n 
niFl'^'SUTK 

TFtOlH .CljU.non TO 73 
'sl.^l^roNfl 3/2 
rUFT=rMli-TfWMn 
fUFTI^ABSCDIFT) 

TF(UTKT.r,T.u.)GU TO 40 

TtlTCH=THTC 

niF'CNsniFi' 

OJ TO 45 
THTCP=T!1TC 
niFTRsDlFT 

TFCorFFl .GT.n.UGU TO 73 

P k 1 T i 5 5 FH I r 

Fi)P4AT(48Xp rnfCKNFSS Ti»» MFTFR= MH1 5 , 8/48X, 18 C IH- ) > 
PHU rVOO, (T(LJ ,L=1,NNJ 
F,J'MAT(IX,RK1S.83 

prtntios 

F ,)P.1AT(2X,128(1H-)1 
RFTJPiV 
F ^ 0 


************$******************************* ****’^*** 


SUBBauTPJE ROSE 


fUMFNSTON Ft?5)' 

C.lM^.n6/Al/F,THTC,K 
FjMMn6/A?/TT,FllN^ ^ 

SlG=-(2.3*1 .?5^1000./TT) 

= »Ffi<4f*2 

FUO=-CThTC'*‘*//l 44)45 

RliTljPti 

Fon 
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C 

c 

r 

r< 


n 

n 


5 

10 

70 


900 

7 

4o 

47 

45 


IN METF!R 


.JORCIOOI 
-TKrrTl,B,ri,rT2,M,TM» 


lOioSfHlATE G/ASc^PLlrLCM^^^CTERlSTTCS HP 
np TPf Pp.'>AT(tPp H^OFfT F AND DETEPHTN aTXQW 

bv SUrjVT.VG time ORPENDENT heat EyHATinM 

OkM) ^hPKiUTiON Tn^niA%ATTn^^^’^ ’ 

’^VAO-’M.TlNr, ti!:mp ir rismdJh^'^ 

n 1 *•<’ (<;’^,s T u "1 T 1 1 C 1 yo 1 , TT 7 f 1 no ■> 

'UVKH.STilN DACSOKCKTSoI vaPftAO’* 

s 

^ ,6.10/ 

DHTi\rb,NA 

4 A'H bO X , n / 1 X , 1 10 n H-i 1 
^A (\ ) = ROO,EMO 
D.) 10, 1 = 2, HA 

D.\ f [IsriAf X-1 ) + i 00 .E-1 0 
OKf 1 1=0.02 
D.) 20 J=?,NB 
CKC JlsCKf J-t 1+0.01 
r.'iofrNnE 
Tl=i09, 

TK,=u , 1 
THTC-45,E-6 
HS=100, 

Rsn=iu, 

AHEA=16.R~6 
PRINT 9UO,AREA 

E'lWMATC'/OX, "THICKNESS IN METER= " , El 2 , 4 , 5X, " M AX 
PRINT 45, TT 

'AMBIENT TEMP IN DEC KEDVIN=' .El 5.8/44X, 20( 1H-) 1 

TC'lND=t 
Tr=,i’T 
CAI.b HOSE 

PE5T=(1 ./STGT)*(THTC/APEA) 

Print 4 7, rest 

FtlPMAT(40X, 'RESISTANCE AT RQOH TEMP. = ' , F8 ,2 ) 
nj 40 T=1,WA 
PRTNTHS.DACn 

EUR AArC45X, 'AVERAGE PARTTCGE DTA IW MRTER= ' , E15 . 8 1 


TIME 


='.E12, 


700 

5b 


CDMpHTATTON OF SIZE DEPEf'iDENT MET.TTNG POINT OF BISMUTH 

ABT=70P.9«"'"'”'’ 

VA=ARl/9.67 

VB=AR1/10, 

Tm=271.+271. 

AD1=nAU)»l,E2/2. 
nC»CVA-(VB*>»=2)/tVB-VA))/ADl 
nH=-f2,l‘54-41/Cl0.2=I^ABT’t'4.184E07) 
rvlB = T«V«*EXt’f 1)B4DC) 

PRINT 700 . TMR . 

FUR 4A1 (40X, *MEL1'INC TEMp, OF BISMUTH DEPENDING ON RADIUS^' .r6.2) 
D!15U U = 1 ,NB 
PR! NTbb , CKf J1 

FOPMATCbOX,'RATin OF PARTICLE SEPERATION TO DIA=',F9.3) 

M=6 

TIME=7.E-5 

Ti?(n=iT 

DU 1 1 Kl=2,11 . ^ 

T12{KI)=TI2CK1-1) 
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1 1 

,r^ 

i.. 

0 

c. 

86 

85 

42 

780 

78? 

781 

788 

770 


01 


00 


Oi,!'" 'ni-tUK ■ 

= ic 

0,)1 still) + 1 

'fl.isu 

'rAF’i>=(TiirC*6.)/Cll + l /CKfjn 
v'a‘(i ) = 7?\pp/io, ^ ^ 

no 86 L=2,10 

no 85' I-=1 ,10 

^«"jT‘^(VAP(I,),CnR(L),I,=l,101 

H 4 ^ i’*i 4 *1“ 1 

TKrM4,E0.1)Gn TO 780 
00 'in ?/() 

SU«! = 0, 

M2 = M 

no 781 K8=1,M2 
K9 = 0 

T4=Tr2(F8) 

KM=K9+1 

T5aT4+( rT2CK8+l )-TI2CK8))*0,l 
TT=CT4fl'5)/2. 

TCnODsl 
CALL ROSE 
T4=r5 

r»uw=suM+f i,/sini)i'o.oi*THic 

TFf K9,LK,y)Gn TO 782 

rOMi’TWllK 

Rl)=SlfM*2,/ARF:A 

Sl=VAPP/(Rn+RS) 

r>V = VARP 

PRINT* nn, TIME 
PRINT 783,5I,SV . 

F :jRMAT( 15X, 'OFF STATE CUPREMT AMP. = * ,E1 5.8, 15X, 'POT. 
lAPRnsS SAMpTiE volts', E15. 8) 

PhTNTUO 
r,U TO 72 
M2=M 
XL=M 

Bl>=XL*0.1*THTC’t'2. 

DMsTHlC-BU 

RMS {RUHM*I.)M) /ARE A 

SU«sO, 

no 90 K8sl,M2 
K9 = 0 

T 1=112(68) 

kusKO+I . 

TS=T4*CTT2CK8tl)-TI2(K8J 5*0,1 

TT=tT4+T55/2, 

TCONHsl 
CAr.L HOSE 
T4sT5 

S|JMsS(jM+(l,/SlGI)f0.01*THIC 
TF(K0.LE.9)G0 TO 91 
CONTI NOE 

ED=fvAPR*(l!+l./CKCJ5)*BD)/(8D»(RM+RS+RD)5 

EMstVAPPl‘KDHM)/(AREA*(R«+RS+RD)) 


DROP 



Q'P 


9 4 

72 

7h 

12 

506 


500 


r> 

K, 


505 


41 

RO 


7 b 1 

50 

100 

40 

n 

•ft 

r*' 


125 


'POT. drop across 

I’KTaTVS^K 

iV?«5?'o!f?6i7Tn‘'So''’ xeter^-.eis.S) 

?Sr'r¥i5i5-^'«>“ ™ 50 

TlMElf}«F;o:hTI«E 

?i:M':'?E'iS^!'{-"”>“ ™ 50 

1.0)00 Tu 41 . 

no iri*‘^06’ 

TU 41 
"" 505 

^ethou^of^rIsectton"' 

TlMt=(TRNtTKP)/2r"””’ 
rifCfrME.GT.TTMXlGO TO 50 
CAT.L TEMEV 

lU 5 05 
M*M»1 

rF(M,GR,2)GU TO 42 
TF(M.E0.1 )G0 TO 770 
fM-sCRUHMKTHXCJ/APKA 
S[sVAPP/fHM4RS) 

Si/ = Gr*RM 
PHTNT*,RM,TlMt; 

PHT<7T 7b1,SI,,SV 

PdRMATt 15X, 'ni\ state CIIRPENT AMP,=:',E15,8,15X. 'POT, ORnp 
IAORHS-S SAMphL VOLTS', E15. 8) 

C'JMTIMUK 
PRI.^Tl 1 0 

E.iPMATCIX,130(lH(=»=)) 

'UMriRiFt: 

s rnr* 

EdO 

StiBHOUTlNE TEMEV 
GAUSS TTERATTOW 

nJ^^ERSTUN^TU nOO)rTT(100),T2(tOO) ,BW(100),TT2C100> 

C.J«*inL/A1/THTC,Bn,DM,lC0RD,R0HM 

r!au,'^!)w/A2/DlF,DlFi,TKIV,TKP,TIME,TK,TTl,E,H,ri2,M,TMB 
CiJ‘4MnL/A3/TT,STGI , , 

PHnC=2,87e03»0,16*4.184E03 
n2s!lThf C*42*RHnC)/fTTMF4TK) 

K5a0 

F 5 s K 5 + 1 

T1 f l)=TXl Cl ) 
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’”H ! t ) = ;rTi(ii ) 
n 4 = 1 ,b74<l/(rit=r7) 

*'*b = U4/7, 

*<'4 = 0 

lOS + 1 

Kj = i,y 
rr=i'ij(K3) 

T c'nls nsS*"' * ^ '* ^ 'f ^ I f i +n (K3 ) ) /2 . 

CAi.ii nnsK 

TX=TTl(2) 

Tjf^.0|^2)TT=(TnC2D+TH2))/2. 

(2)+D5»Tll C3) + CTHIC*E)^*2*SIGT 

T‘UnTfiof 

r af.r Ans ) TT= c tt i f i o ] +t 1 ( 1 o ) ) /2 , 

'■^iMriO) = U5*Tl1 C9) + Cn2-D4)*TTlClO)+D4*Tn (11) + CTHIC*E)**2*SIGI 
n'KlOJ = BMC10)/(D2+n41 
Rt=l)S/CL)2fD41 
TFCK'^.GF 2)Gn TO 205 

Tir2) = Ol\lC 2 ) 

D.j ao K2=i,if> 

Tl rK2) = BMtK2)+Rl*T1 (K2-1) 

00 ru''n'rNni.; 

KB=0 

205 KH = KRn 

nj qs K4=2.io 

TE('K4,r,li;.10jGU TO ooo 

TFCK4.l.E.2)Gn TO B05 

T2(K43 = HfUK4)+Rl»(Tl(K4-n+fl(K4 + l)) 

GO TO 05 

005 T2(K4 )=6N(K43+R1*T1 (K4+t) 

GO TO 05 

000 T2CK4)=hMCK4)+R1*T1 (K4-1) 

05 noMriNnc; 

T2fl)=Tl(l) 

T2(11) = n(tl) 

0.3 lOy K4=2,10^ ^ ^ 

TCHE»AnSCT2(K4)-Tlf K43) 

TFCTCHF:.GT.O,l)Gn TU 108 
10O COGTTMIK 

GO TO 110 

10R n;j HI K4 = 1,U 

TlfK4)=r2tK43 
111 cuMnwoE 

G.i TO 105 

1 10 00 120 K4=1 ,11 

Til CK41=T2(K43 

120 continue „ 

TKf K5.I.E.y)Gn TO 125 

0 IFsTl CWl-TMR 

Of.i’l = AH5CuTr)^^^_ ... 

TFfOTEI.I/E.l.OlGO TO 406 
TFCuTK.IiT.O.lGn TO 400 
TKPsTlME 

GO TO 405 
non TKMsTJME 

no5 no 12 T=1.11 

TIU T) = TT2CI) 

12 COMTTNOfc: 

GO TO 4U7 

406 nu H T=1 ,ll 
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1 4 

1 

4. 07 


12 


T], 9 ( T ) r'i 2 t T ) 

Cl “I I r oHL, 

f’KTtMT US, (T7(K4),K4 s:1 1 i ■» 

'nsTunf; ' ‘ I 

F.'jrv 

Kic-M-'n jjT (, Wf.,; priKF 

COMDUCTiViTi 


TFUCIUU.^O.IUGU Tn 72 

RKTUVif ^ iU.)/SlGl+D!«/ ( 1 . /ROHM) J 
Sr-^'n 
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r 

ft 


15 


70 


1 io 


80 

90 


70 


50 


c 

60 


1b 

100 


10 

50 


^'U’0iAT^lJ^l np 


C2^ TWO PROXIMATE. 


sriKpr^srsiiu. 


Ot-- A« Pl.RCrkTC FTPT.D 

* ii( 4« A *4 J * 4. Jk .t a, .1 ». T 


AXi5,B=sEMi minor axis 


HAIHI 5) ,CKf 15) 


R.'UH n = 100 .E -10 
no 15 KKs2.N1 

rArHKts)=RAf)(KK-i )+50.E-10 
rOMTlMIE w.n. iu 

CK(i)=0,oi 

nj 20 0L=2.N2 

^^(LU=CKCf,L-l)f0.01 

riNirhUE 

fV,) iO KK=1,K1 

PhTnT r20,RADCKK) 

FjPMrf(25X, 'RADIUS IN METER%E15.8) 

i 

RRTrrr 1iO,CKCLL) , 

P'jf'-jM'f 'SEP. TO DIA. RATI0=%F5,2) 

'n‘''«E=! rME+0,1 
CJHTIMIfc 
OAF. I, SWDYN 

IFfSMiN.I/r.O.IGO TO 80 
TPCADIFI .LT.0,1 )GU TO 45 
Tf (AOIF.DI' 6,0)GO TU 50 
TlME=TrMF-0,{4TIME 
OAI.O SWOYn 

TffSMIN.I.T, 0.0)00 TO- 80 
Tr'f ADlFl,l.t.()-1 )G0 TO 45 
TFf ADIF.T,T,O,0)Ga TO 60 
Go TO 70 

TIMEsTTME+O.HTIME 
CAF.O SWDYN 

rFfSMlN.ia'.0,0)Gn TO. 80 

TFC ADlFi .LT.O.l )GU TO 45 
FFf MUF,GT.O,0)GO TU 60 
GO ro 50 

METHOU UF BISECTION 
T1ME=( rKN+TKP)/2, 

"’Ar,L SWDYN 

TFC ADlFl .LT.O.l )G0 TO 45 

GO TO 60 

CONTIMFfi 

print too, (TTMF.SKAJ.SMIN) 

FJP^ATCbX , ^switching TTMFs 
l£t5.H/2X,127(tH-)) 

C OMIT ONE 
CONTTNlIb 
SrOF 
F MO 


,E15.8,5X,'R=',F6.2,5X, '8= 


I 




■?:^f;S?'iSrSr“SK«cs of s.ttchino 


DEMENSION 1 ”l 9 1^PHTF{-1 •12.»t*12)!,U(*’l!l2,"lSl23 

niMh'JSTON ^ I 

!) l''^E?>ISTUN IFX(-1 :15) , AUYt-i .p) » YC 1.12, 1.12) 

COM tlON/ A1 /cot ,C((7 ,PHII , Al ,B1 ,FUN 
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1 0 
lyn 


■?o 


b 

1 / 
15 


25 

?(> 


40 


' - I . +IK t I,,Ij) 

:'! = 4 , ♦'S'fTKi'j/ CPPt^KPU^RADCKK) 1 
K = .50'Vr ( C * t P-SORT f P) ) ) 

!’y| = f .ElO 

PbrpX(nj=u,;PHlBYCO)=:0. 

iu c p ) = 0 . ,* V xp ) = 0 , ; nx ( 0 ) = 0 . j vy { 0 ) = 0 , ? XX ( 0 ) = 0 , 

*>ij I J 1 = 1 # 1 0 

xxf n=:xx( c-nto .1 

oiTBxtr)=pmBX(i-n 

pntiV(n=pifiB)£(i-i) 

’jxn i:!iuxa-n?vxa)=vx(i-i) 

;] f 1 } = ; J J f i - 1 ) ? V Y c n = V y ( I - 1 ) 

f,!so 

Pn'n = kAU(KK)*l,EOB 
\ .f*;o0 

Wl = t4*J,K;0B 

M,X10J=I.U 

r so * xcf) 

S0 = CKCt»Li)*2.*RAul 
5 = .50-2.-t‘CA3-RADl) 

KsE+f ftO/B) 

C l = t0.75*EP*Epn*F**2'>/nE?lS 
r2 = S!IlRM/0F:MS 

C i=:2.i0i»ci *ALnG1 0(3./(3.-X»’<‘2) ) 
r4st Ai/»3j**2+ri.-X**2) 

C5*l AJ/H3)*(C4 + 1 jfX>t‘=t'2) 
rb=:(r4tx**2)+*i .5 
r / = (.C5/C5)*(1 ,/B3) 

Cd=(C2*l.K0B)»(C7-2./kA0l ) 

C9=,54f lIXtT)t*2+VX(lU*23 
rinsC3-cfi~cy 

PHTBXTf .l) = PHTBXtT)+0,14fTME^ClO 

Y=.5oP’i(ls-XXCi)’t'=*=2> 

! Ff i jEO.njGu TH 35 
Di 36 3=0,10 
rKV.GT,XX(J)K’0 TO 36 
Tl=J"l , 

Gi,) T'O 37 
CJMTTivfnt’ 

AOKtTJ=CY-XX(It))/H 
ra'^TTNOE 
T = T 1 1 
x=x+o,i 

pl'l f a‘/j)=PHT6XTr5) 

cafn’Tf'.ufe; 

rONTIwUE 

A6Y(»J=1.0 


cS=C/i3/rt3i*(Mp>X*»2’ 

i! li ; Sfe....... 

rys:5icU¥n?*»2+VYfI)=!'*2) 

riO=C3-C8-C9 


adif 



'r‘J V i ^'1 a ^ 1 * f 1 0 

f ^ C X # t’.j> 3 m 0 } Gij fO 46 

^■'J 4S ,T=n,J 0 
TffX.GT.XXfJU GO TO 45 
T I =:J“i 
TO 47 
CiJM'j XMJK 

AhYt T)=:{X-KX(iT))/AK 

rilMTTrjnK 

TsT + 1 

y-YtCM 

Ir'(.l,LK.10)Gn TO 40 

r-’YALOATlON UF l.APtACF FQUATION BY GAIfSS Ti'EKAl'IUM 

®Ste Sl^ ’RSW iSA 408^ 8SB 99> gggf |Wi ifgg^ 9ltS 0811 900 1801 0118' m0> 090) 40! I0K 4988 MRK 9 

Ai=l ( a/bX*H)**2 

iU = AKtt2 

C 11=(l./Al)/(2.*f l./Ai + l./RU) 
r:o2=(i./Rn/(2,»( i./fil+u/BU) 

Li = 0 

N4=9;wB=4;K=0 

r.l=Lif 1 

1)0 50 J = 0,0 
K-K + 1 

rFCK.GT.5jG0 TO 55 
D^=Ai.nJ);ox=u 

TFC J . LO.Mb)DX=ALX(N45 
00 bO T=:0,M4 
CAr.L FTIHFAY 

OriTK J , J)=FUN 
r.joi'TfjnK 
GO in 5U 
M 4 = H • fl 6 = 5 
rFCK.GT.bjr.D TO 65 
n*=:ArjY(J);DX=AliXCN4) 

DU 70 1=0,54 
CABO FONFAV 
PHT F( 1 , J)=FUN 
OiJOTTWOF 
GO TO 50 
54 = 7 7 M6=7 

TFCK, GT,A)Ga TO 75 
n5{ = Al.i( J)?DX=1.0 
TFf J,F0.Nfe)DX=AbXtN4'> 

00 BO T=0,M4 
CAGU FUnF.kV 
PrttFCI, J)=ru5 
COOriNlIE 
GO TO 50 
'*!4 = 5?56 = P 

TFfK.GT.PjGO TO 92_ 
nx = ALX(a)|DX=Ar.X(W4) 

00 b5 1=0,54 
GATO. FOOFAV 
PriTFf I»0)=FU5 
COMTTj'jiJF 
Gu TO 50 
y4,=:4;i'i6=9 
no yo T=o,H4 
nx=ALX( l),*nx=i. 

TFCI ,F0,N4)DY=ALYCN6) 

Cwr.b FU5EAV 
PtiTFCI »J)=FUN 

coyrronE 
cowttnUg; 
no 95 0=0,9 



1 00 
1 OS 


120 
1 10 


1 30 


oin 


620 

600 

n 


1 /o 

tb5 


IhO 
1 76 


190 
1 8S 


li-’f J.i%0.«JM4=5 
OJ 1.0u I=0,,N4 

CJOir,',ot; 
oj tn 130 
r::}'>ITr^t?K 
o.j no j=o,9 
0 4 = 4 

[[;:[ j*;i5«5{;;|=^?IP'('T.t’0.5)N4=8;TF(J.K0.6.0H.J,E0.7)W4=7 

n.j 12 u'I=6,S4 
Oi:ir|a,j)=PHrF(i,ji 

V ’. j j t f fli U tl, 

CUOTT.VOe 
aa TO 140 
CiH'TtfS'IF; 

Hij faOl! J=0,q 
fj'ls'l 

= J.fe;0.6.0R,J.E0.7)M4=7 

If CJ.fcg.s J04=5 

O'J 610 i = 0,W4 
PHTIC1,J)=PhIF(I,J) 
cj^rrifMoe 
ao TO 6U0 

620,CphI|(T,J),I=0,N4) 

F'jOKAl'dOtU.S) 

roMiTNUL 


CJ’^PUTATIOM OF velocity COMPONENTS 

Mi=:yyH6 = 4;K = o" 

00 160 J=l,9 
K = Ktl 

rFCK.GT.NblGO TU 165 
ni=Ar.nj);nA=i.o 
rF(J,t.0.Nb)DX = ALXCN4) 

00 170 1=1, N4 , 

CAIO. VELEAV 
'1(1 ,,JJ = SUM1 
Ytl , J)=S!!M? 

CONTTNiii^ 

GO TO 160 
M|r8?i\i6 = 5 

TFCK.GT.NoIGQ TU 175 
O^sAI.i ( J)J0 X=ALa(N4) 

00 190 1 = 1, Ai4 
CAr.L vFLEAV 
lu r ,jj=6'JMt 
Y(r,.T) = SUM2 
nOMTIOOE 
Go rn 160 

0..|:=7;i,6 = 7 

TFCK.GT.NDGO TO 185 

r)V = AT.lcCJ);DX=1.0 

TFCJ.tQ.Nt))UX=ALX(N4) 

no 190 1=1, N4 

CAt.L VELEaV 

n(r,oi=snHi 

’/(T , J)=S!J«2 

CONTTNOE 

nj TO 160 

W4=5?N6=8 ^ 

TF(K.GT.W6)G0 TO 195 
Di(=Al-X(JI,*0X=ALXCN4) 

DU 2001 = 1 ,M4 


21 6 

n(r,.i) = snMt 
'nT,>IJ=SUi« 

700 rai'^TrfMiifc. 

'Tu 'I'l'i 1 60 

lyS :-! 4 = 4 ; 1.5 = 9 

''I') 210 i = l,iv4 

9a = a?.x( i) ;nx = i , 
rt*’Cl.ii.O.N-'nriY=ALY(N6) 

MKlF.hV 
vt I ,.TJ=S1JH2 

?io c j'‘i n'Ktii:; 

160 C-JOIlNOb 

no 710 J=i,y 
N4 = 4 

GO m 710 
710 GOmHm‘!£ 

^ ™ o"¥”nfI"'BnUMDAHY” 

'\2 = A*!i 
B2=B*AK 

n V i 2 u " J“ 1 ^ y ^ ^ 1 ^ ) -Pf^ If’ C 0 , 9 j 1 /b 2 I y Y 1 1 0 5 =0 . 

Tb( J.b0'.6.6R,wI,EO,R)G0 TO 230 
lKCJ,GT.4lGn TO 400 

jY.u+Yai-soRTci.-cFLOATCJi^o.ni'^a) 

0|.S2=( AOYC J-n*H-DTSn 
TKC J.GT,0,0.ANn.J.T.E.4)i44=9 

PHr = nJTS1 /f ALTCJ"n»Hl)+PHTFCN4, J-1 )+PHIBYT(J-t )*DIS2/CA1.YC J 
t ) H ) 

OVf J) = CPHIHYT(J)-PHIF(N4, J))/(A2*At.Y(a)) 

Vlf ( J) = (HH1BYTCJ)-PHI)/B2 
GO TO 220 
230 CaMTToIIE 

IFrj,h0.6)no TO 240 
■/Y( J) = (FHJB^T(J)-PHIF(J-2.J-n)/B2 
Of ( J) = (PHIRYTt.J)-PHlF(J-3,.T))/A2 
GO TO 220 

240 V<f J) = (PHi:HyT(J)-PHlF(0+7,J-U)/B2 

Of ( jl = (PHlRyTCJ)-PHIF(J+l,J))/A2 
Gl) TO 220 
400 GOOTFrOE 

TFrj.E0.7)04=7;lFCJ.FQ.9)N4=4;IF(J,EQ,5)N4*8 
nisi = FI,0 AT ( M 4 + 1. J ♦H-SORT ( 1 . - ( FLO AT ( J ) * AK ) *»2 ) 

Pofi(UTSl/Oj*PHIF(M4,0-l) + rDTS2/H)*PHIF(N44-l ,J-1) 

0»:(J) = (P0I8iTtJ)-PHlF(N4, J))/A2 
Vf C0)=(PHIRYT(J)-PHI)/B2 
220 rOOTlwO'E 

VX(0)=(PHXfiXT(0)-PHIF(0,9))/B2;UX(0)=0. 

OX(l0j=tPHTBXT(l0)-PHiF(9,0))/A2yVXCt0)=0. 

no 250 T=1 ,9 

TFC I .E0,6,nH,UE0.8)G0 TO 260 

n l S 1 = SOAt 1 1 CFLOAT (T-13=*'0.1)**2)-SQRTC1.-C FLOAT C 1 3 *0 , 1 } ♦ »2 ) 
'nB2=:ALXf I-ni'AE-DTSi ^ „ 

P0l4?Ul5l/(ALX{r-i)lAK)3»PHIF(I-l,N4) + (DIS2/(ALXCI.l)W)3* 

1 P H f X T 1 1 •*' 1 J 

OXfn = (POlBXTCX3-PH|)/A2 

VX( I )s(PHIBXT(I3-PHIPC I»N4>3/CAI,X(T3»B23 

GO TO 250 
260 COOTTnUE 
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2/0 

'130 


250 

r » 


R ()0 


QOO 

350 


Ti^’U TO 270 

^ - f JJ » ^ ( n -PHIF ( T , T + 1 ) ) /B2 

'J !( f J 1 = ( PIU RXT( T )-PH1;F(T-1,T-.2))/A? 
^ ^ J. ^ t ^ X T U J - F H I F ( T , I - 3 ) ) / B 2 ^ 

F j^TT - artfc ; 

n 52 -AK ni ^ i ^ ^ 


1 r : f P W I R ¥ V f T 1 e-Dti '¥ t t i ? 4 h / 1 % n 


c\h i i H Ue’^^ i a ) -PHirt 1 ^W4 ) ) /B2 


r ;jWFU'i’A|'tuN HF NKw boundary 

>VF = Hf AK=fTlME*VX(0) 

A = RAlHf<K) + f.HAD(KK)/RF)*^2 

IVxHf 

nj BOO i=uao 

OHTbK(.T) = PHIBXTCT) 

OuT « aCIJ = PHlB 7 TiT ) 

rjNrTwut:: 

.S''.AJ=:A*1,K5 0/SWIN=R*1 .FIO 
TFCB.l.T.O.OJGU TO 350 
TFCL.LT.KDGO TU 300 
ADrF=;iMA0-S01 
ADTFIsAHRCADtF) 

TFCAOlFl.LT.O.DGU TO 350 
Te‘(AniF,LT.0,0)Gn TU 900 
tko=txmg 

GO TO 350 
TKN=TlMfc; 

F^fiTURN 

F^n ! 

^.^ iifj ^ fi . i ($ it : ilf :: iisi . i ^ ilf .^ 3 ^* i ^ t ******$'******^*******************^****^*’**; 

SURKnUTlNK FUNEAV 
Fi/AUFIATION UF PHI FUNCTIONS 


OIMKMMOM PHTl("1:12,-l:12),PHTBYT(-l!l53,PHIBXT(-l:15) ; 

CQMf.UliWAl/CU1.,C02,PHli,Al rRifFUM 
rO’Hinjsj/AI/I, J,PHlbYT/PHlRXT,M4,N6,DX,DY 

iF ( tJ - ua . T . ojPHTifi , j - n = PHTia , j + i ) I 

TFfC r-l) ,UT.n)PriTi(l-l,J)=PHIl(i+l ,J) 

Tf f J.FOjiGO TO 15 

Tpfi pn W 4 i('^n TO IS * ■ 

FuH=CUl*rPHftCT-J ,J)fPHII(I-t-1 .J))+C02»(PHlI(T,J-l)+PHini,J + l; 
hii ' I'O 1 c ) ' 

15 n = 7pl t ./tAl*UY) + l ,/(Bl*nX) ) 

F4i;f2,/CAl»(DY+i .DI/O , 

CA2 = (2./(/i1»9X*CGYtK*U^^“ 
r n - tx./oii + cux + i . ip/o 
ry 2 =( 2 ./(Ht»nx^(Dx+i.)))/u 
lFf,J.K0,0)G0 TO 20 
GiJ 1.0 3U 

20 ‘’ iin ( l,J + l )* PHiaXT ( I ) 

FFf t *K0.N4)PHn(Tt1 , J)=PH1BYT(J) 

Gtj rn 25 

30 PHT.ICN4 + 1,J)=PHIMT(J) 

1 J+134CY2 

10 RKTURFI 





«,!n|.:ntn’i.M[-: MW.hW.^v 

?!v’Ai,"irnntt OF velocity COMPONEMTS u,v 


n I '^Kf',sTLiM PHTFC-1 J 12,-1 ;12),PHiaYr(-l;15) ,PHTBxrf-i : 15} 
"'.!•< 1/T ,.T,PHTByT,FHIBXT,M4,i<6,nx,DY 
r t‘ihnri/A2/A ,H, A^,H,PHIF,SUMi ,SO«? 

FF1-1 ,/t2.tAl-ii) 

CF2=:t ,/C2,tAF + B) - 
TFf J.EO.^ljno TO 15 
TFCi .t',0.fJ4)Gn TU 15 

no THIO 

r:Ai=i ./t!\=iHi:nY*(nY+i.)) 

CX'^rsnY/tAl'Hl^CDY+l ,) ) 

CX'l = (U-Oy)/(A*Hl'UY) 

CYlsI ./tB»AK*l)X + fUy+1 ,}) 

C:t'2 = nA/CF^=«cAK*(.DX+l.}} 
r‘(.i=(i,-nx}/(Hi=AK^=nx} 

UT J.h.Q.9}Gu TO 20 
Of] TO 3 0 

PHTF(1,Jfl)=PHTBXT(l} 

U'C 1 .L0.N'4)PH1FCT + 1 , J)”PPIBYTt J) 

Gu to 25 

PHTKf i.‘1 + 1 , J)=FH1F^YT(J)^ , ^ 

U'f J.E0.M0 }PHIFiW4,J+1}=PHTBXTCn4) ,, 

piiMl=CX3 *PHiF(T + l .‘JXX'i^PHlFa-l ,J3-CX|*PHIF(T,J) 
SUM2=CYI*PHI.FCT , J+ 1 ) -CY2»PHIF (T , J-l) -C Y3*PH1F ( I , J) 
RETURh 
KNP 
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C 


650 


6 

5 


1 b 
1 70 

1 10 


PE05i<?.^yE.TO_Fi^0 iHe.optTKnM^ 

__ _ , .. ,_ ‘thpory 

m\r SINGLE STRRnO CHAIN 

Opt. t y J =;uP.pni,AKT6ATTUP TRTpT.Ef FOR SPHERE 

glass Sairxx 

iJJpfi'SijXHUuKXLaoiriocroKXHaoirmlanxccro) 

R P no), Yt 10,10 ),ARCOEFC100).DFP(b,b'),IXALPC 100) ^ 

XX(300),ALHD(i6o),FPBl(100),EPB2ClO(n,AliCOFC100) 
CSiJM(lO) ,CPEFUN(10) 

N^*W4,XX,ALMD,EPB1 ,FPR2,0EP,EXALP,X,I,ALPHA,snHt ,PEFU 
C j-l.n'Jiv KK, AHCuF,AbCOFF, irfjND,<^Pd, rB,VFP,WpX,FlF, AVD.FPMEn 

M-2 

Nlxw+l 

02 = 0 

WpB=(l6,n^l.602E-l?)/l,054F-27 
Tb=l .064R-27/ 12,8 38 1»1,602F-1 2) 

VKH=,l.822E08 

4px=4[>n 

FlK=.y4i 

nn^'A A,H,C/l ,(),0,5,2,0/ 

>11=1 ’),*N3 = 189 : 

nEAI)’(‘, ( ALMntJJ) ,FPRI (JJ) ,EP62(JJ'» , JJ=1 ,N4) 

Pnf^J T^, (ALMDCJJ) ,f;Ptil(JJ) ,FPB2(JJ) ,Jv1=I ,N4) 

REAUl‘,{XX(TI),TI = l,iO) 

PEAU», CtnEPtM,Mn,M=l,3),Ml=l,2) 

PKT,aT%, CCOFP(M,Mt) ,M=l ,3) .»-U = l ,2) 

READt, CtYtJ,I) ,J=1,M),T=:1,N1) 

READ*, (EXAJ.P(JJ) ,JJ=1 ,N4) 

DJ oSO JJ=:t,H4 ; 

ZAPP=KP81 (JJ3 , ^ 

FPRUJJJ=EPB1 (JJ3**2-EPB2(JJ)**2 

FPrt2(0.T)=2.4EPB2(dJ)nApP i 

CiJ.MrTLHE 

PRTwT»,CXX(1I),1I=1,N3) 
nj 5 I=l,wt 

0',) b ! 

XlT,J)=Y(U,n ; 

PHTH r)t‘,X( l,J) 

CuPTIi'itTE J 

cjNXTi'^oE ; 

STMPLEA UPTIMTSATION lECHNlOUE 1 

K = 0,» T.C{tND=0 ^ 

T=1;RK=.l 

CALL EHNEAV 

RA=0.3*SlJMl/(-l.^PEFnN) 

rcnun=Tcni^n+t . . 

OJ 15 T=1 ,N1 
CATiL FUNFAV 
FITJsAI.PHA 
CS'tMd 1=SUM1 
CPFFHNf I)=PEFUN 
CONTTMIE 
K = f( + 1 

xyso 

K y ~ K 9 + 1 

F<AX=FC1) 

TF(tSAx!GLF(I))GO TO 20 

F.4AX=Ftl) 
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7 0 
72 


75 

26 


4 0 


30 

15 


45 


60 

50 

55 


65 


75 
7 0 


QO 

90 


410 


! )-J 


27 

XafJ)=X(TH,J) 
CU^^ilTnni;, 
FiTH=i< (1) 

TL=:1 


fVl 25 T = 2,M1 

•JKfK'»ilO,LT,K n)GO TO 25 
f’ h I .iSf f I. ) 

c.)*' Tiiuit: 

O.J 26 

4 i=r.i 


ou 30 J=1,M 

.5!)^' = 0, . 

no 40 T=t,Ml 

Iff I .KO. THIGH TU 40 

S)J-’^ = Su«tX(I,J) 

CijWT TfillL 
HtlWsPUfJ'/ftN 
Xtif j)=SilM 
r J^iT^wHfc 
n;i 35 ^J=1,M 
X(n2,j)=X0Cj) 

T = M2 

CMf.l. HIHFrtV 
FXOsflGPria 
SUMOsSHHI 
oa 45 0=1, N 

Xuf J)=(i,+A)*XOCJ)-A*XH(J) 
X(02,U)=XH{U) 
r JNTTUHK 
r=M2 

nAf,L FONFAV 
FXP=Ai,PHA 

TFfFXK,LT,FMTN)GO TO 50 
60 1 = 1,!^! 

rFd.hO. THIGH TO 60 
1F(KXk.LF.F(T))Gn TO 65 


FOOT fwOK 

Gu TO as 

CHNTTMPF 


HO 55 J=1,N 

XE(J)=C*XK(J) + fl,-C)*XOUT) 
XtM2,Hl=XECH) 


CAr.L FUHFAV 
F>CF=ALPHA 

Tf CFXk.,LT.FHTN)GO TO 

CJTVTMHE 

H!! 75 J=1,N 

Xlif J)=XRf J) 

x( rri, j)=XHCj) 

G'J TO 8 0 

C.J^'TlNnE 

nj 90 J=1,N 

X3fJ) = )(EfJ) 

XCTH,U1=XH(J) 

C-l'iTTwdK 

rv-| 95 T = 1,N1 

rKd,fc.9.Iii)G0 TO 410 

Cn) TO 95 

Car.o FONKASr 

cs0M(i)=snMi 

rpKH3H(I)=PEFUW 

F(r)=AOPHA 


70 
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Q5 
1 45 

1 1)0 


non 

700 


85 


120 

1 15 


125 


1 ,35 

400 

140 

130 


155 

150 


750 

105 

1.50 

4,30 

190 

900 


.S'i'-ist), 

0. 1 10(,i 1=1, fvl 

F(i •=r>ii'‘i+ fFC n-Fxo)»*? 

"0 ;y=^0KnsuM/Nii 

1. 'u hO<* 

v'iM)=xrt,j) 

" J i'l TivUr, 

If (kO. <71^150)00 '^0 500 

" j‘' 1 

ffc-'f FXn,ijT,FM4X)Gn to 115 
OJ 170 J=l.H 
X5.(J)=X5fJ) 

'v.) 125 J = l,Ki 

^ 1. ''2 
T = 2 

FATjI, lOwFAV 
FXC=ALPHA 

TK(FXC,GT,FMAX)Gn TO 130 

nj 135 

xhcj)=xc(j5 

X( Th, J)=Xrt(vJ) 

no HO 1 = 1, N1 

TF(1 .t-JO.TrilGn TO 400 

G.l ffl 140 

7Af.L rONFAV 

r.siiMf I )=suMi 

roKH'HrKi )=Pb:FUN 

Ft T J=AI,Pf1A 

CuNrJi-jiJF 

C,-J in 145 

CjoTiritre 

no 150 1=1, Nl 

no 155 J= 1 ,N , 

Xtt , Tl=(XtT,J)+X5tJ))/2. 

CJOTINOF 
rjMi'iNOt; 
nj 750 1=1, M 
CAt.L bllyF.AV 
Ft T )=ALFHA 
GS'lMt n=SUMl 
CFFFHAtX)=OEFUN 
CJOTIMIE 
nj TO 145 
Cj>''!'rTM5e 

TFCK.FO.I )GD rn 1«0 
PiJ‘^2 = SUM0 
GJ ro 190 
FU'n = SIIr>10 
OXSOHI'HK 

n,j 4 30 ,1 = 1 , Nl ■ ■ 

Ft f HCSUMtTl-RKl'CPEFnNtl) 

CON,l’IN!JE 
no TO 170 

nr.j7=ABSt (F0N2-FUN1)/FUI125 

r'i3rNT + ,UE7,RK 
prints, tXOCj),J=l,N), 

prInt'^^OO* ( ATiMDC JJ) HBCOEFf JJ) »EXAIiP( JJ) , ABCOEtJ*!) , JJ=1 
FjR,^iA‘CC2X,4(4X,F10.5)) 

P ^=H K ^0 , 2 
PkIPT 050 


,N4) 



nsn 

1 4 0 

70 0 
f>UO 

?b0 

BOO 




1 0 
15 


10 

70 


35 


2«3 


sKOri? 

nj 4-10 l„ = l,oi : 

f ( U=t:3!lM( T j-HK^^CPEFONf 1) 

! ;,i!t 

Gj rn I/O 

fi.iT OUn,FlF,A\/D 

n ' t n I ^ ^ V f 'r % ^ ^ ^ ^ ^ T ^ S I ZE= ' , F 6 , 1 5 

Pnlr^ I ^ ^ i XU f J) , J=t , i ) 

r‘il :: I . r I' Y ' i ^ Z- ? *” ( J J ) , ABCOE ( J J J , J J= 1 ,04) 

> i P f'i A i C / X r ^ t 4 X # H. U « 5 ) } 
r^rnp 
K , n 

f 4^4 ^ 4 tils 

SiiBHriUTIMF. HmEAV 

OF OPTICAL ABSDPPTinM cngFFTCIEiT BY MOPVS THEOi 

ni MF'fSTpO Xxpoo) ,AbMDf 100 ) ,EPhI7Iou )*EPB 2 (100)7^(100) 
nr K nsTD jfi t onj,H tioo3,w(ioo),Fpsicioo),FPS2(ioo) 

ni''f.MorilM PEPaUI 00),PFPS2(100) 
ni '^K ihTllN PEP(5,5) ,|}PL(b),ALPHU5),ALPH2(5) 



-, 1 ) 

rvl=X( JL ,2) 


AXI=(1.-{ZT+SI)) 

TFdPooo.eo.ncio m ro 

Fpytn=i.o.o 

Fpo] =f,PMFD,*KPi^2=0. 

OPT.U 3 = 1 ,/3. 
nj 5 M=?,i 
npr,(.M j=uPL(M-i) 

AvA=a, 

30*1 = 0. 

B-4AX = XX(U 

n.j 10 ri=2,M 

rF(XX(Ii) ,CT.BMAX3BKAX=XX(T1) 

C')MT1N(JF 

ffzO 

A = 30, 

HzAtS. 

K=K + 1 
A!F = U 

03 20 13=1, {v3 , ™ 

fF(XXtTX) .GT.A.AND.XXtll) .LE.B3G0 TO iO 
GJ TO ?U 

'JK=i4F + 1 

OMTThUE . 

Xl (K)=(A+a)/2.*l.E^8 

’.tpi (K3 = rJF; ANF=f1F 

i C K) :sj F 1 t K ) 

Ai/A = AVA + FLCIAT(NFUK3)*(A+8)/2. 

Si,P'i=SUM + Xi(K3^*3*ELUAT(N1' 1(K)3 

rF(A,GT,BMAX)GP TO 35 ; 

c ' 1 T n l S ^ 

Ai/A=CAVA/FbOAT(N3))*t .E-8 j 

A»/I)=AVAM ,e08 I 

3f T) = (?~if2L/7.)*3,F103/(ALMD(J)1'l .E-43 i 

Gl:stCAVA+2.^TB*VFB3/(TB^AVA3)**2+!«(J)**2 . 

iH = w( J)>*2 + 1,/TB* + 2 ; 

FP32(^4=FP82CJ)t(SpX/2cj))?C«PX*lQRl(Jn/Gl)-(WPB/(TB*WCJ)> 

i‘*(;^pa/Hi 3 ; 
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Qij 


4 6 


50 

45 


50 

55 


70 


R5 


75 

C 

■fW*' 


TO 48 

^*T6»VFB)/(TB*X1 (L) ))4»2+W( J)4=»2 


■(wPB/(TB*ri(J) 


50 M=1 , 3 

-‘‘' oOi til* P l-rKPSl cr, ) *DPL (M ) 


C'j ■' i'fwtfh: 

r^l'ITT^OK 

'Ol J = l,^i4 

^l=u 
A ) r i ) , 
r< I ^u. 

Ki rK1 +i 

>'‘4 45 I..= 1,K 
^c'Oh I .Oh, 2) on 
■41=qxi if.Jf2v 

:>1 -=0! - a( JJ + ^7 

pf;p;--.i (r.)= . „ . 

1 f .V P I i / 1 1 1 f ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ° ^ ^ ^ ^ ^ ^ 1 ’ ■ 

4hpi'-0^^^ U.)**3=!'F1CL))/SPM 

4!/’J=ur 

^p=p|^;pfu (5)«ppMi 
f u)-RpM2 
na 5^' M = 1,3 

r> p - ■' - - - 

np: 

A 2=:CPt'!»'2 + UP+ *2 

‘ ‘ 1^’ p \ U . / p p c AP *CP+dP *DP ) / ft 2 
A(.P2'=ALP2t { 1 ,/3, )¥fBP4rp-AP + Dp')/A2 
Ca''M' f iv'ir.; 

41 "/\1 ♦AH1.!)4A5PJ 
ni =ia jCl,.)4ALP2 
Ci.ifi n wn>,; 

4l=B’P5l (JJ-EPMI 
Bj =i',:Ph?caj-fc;pH2 
nj 55 Hi=l,2 
A5PnUMn=n, 

AbPhP (Mi 1=0, 
nu oo 4=1,3 

r i=ii;PMl *a,-DKPC^,MU) + EPSt CJ)*0FFfM,M1 ) 
ni=EP'*2«'(l,-DEP(M,Mn)+EPS2(J)=t'DEP(M,MU 
A4=C 1 ♦*2 + UT**2 

AI,.PHUMn = Al.PHl (Ml) + f l./3.)*fAT*Ci+aT*DI)/A4 
AUPHPtMlIsALPri^CMD + tBI^ri-AT+DU/AA 

Pir'l’TOtJE 

cjin'Ti.at 

R I = 2. 1 * A 1 + S 1 4 A bPH U 1 ) + AK I * AT.PH 1(2) 

= KSr^ALPti?!! )+AKI»Ar,PH2C2) 

A = t ,-(■(.?. /3,3tKIF*Ht 
B=(2./1.)*FiF4R2 
r = 1 ."Cl ./i,)*KTF*B1 
n=(l ./3.)*FIFtR2 

KPMKi IJ)=E'PMFU*1(A*C-B’»'D1/CCT^2+P*»21) 
rpMF.PC J)=:KPMFi)*((B*C + A*Dl/CC**2+n**21) 

ac;VJ=AB5fLPnl-FpMEl tvlj) 

l)EV2 = AHS{r:Pi42-EPME2CvJ)) „ ^ 

TFfDFVl .l.T, 0.01 .amp. nKV2,LT,0.01lGD TO 70 

Fpwi = i-;pf>iFi(ji 

=:<»'''2=BPMF2(wn 

^ ; I ''I* I I C)! 

A.>Fi)FC4)"(W(JJ/3.E10)*fEPME2(J)/SQRTCEPMFl CJ))3»1 ,F-4 
AoCLiFFr J) = ABCOE(J) 

\rtrLSEC‘T)=:AT.UGlOCAbCOE(J)) r .v r,. , .v 

I F(h;f»fiF;i (J),!.T.O.)PRTNT*,Al.MO(J),EPrtEl(J) 

Cj'flTOllE 

SU'-UsSOmI + (|xALP(J)-ABCOEFCJ))*»2 

FaPMOhAXTUM OF INTERNAI. y”?Tf9f T'pf-V^l 

PF;Ftifi=-l ./XCI , U-1 ./X(T ,2 ) + l ./CXCI , U+Xl T » 2}-l * ) 



/» ! , P ; i ft = S lift 1 - R K »PEF UN 



Addendiua to section 2.2.4, chopter 2. 


In section 2,2.4, 'it Is true that we did not 
discuss the ' justification of the boundary conditions. Only we hm'c 
referred to T^, Kaplan and D, Adler’s paper ^Hef, O'?, page 19? J 
where similar boundary condition has been used for chalcogenide glnrsy 

material sandwiched hetwccn two metallic elect The Justification 

of this boundary condition is based on the an sumpti on that the 
metallic electrodes have much higher thermal conductivity compared to 
glassy material so *that the metallic electrodes can be treated at 
infinite heat sinks which fix the end temperatures at the surfaces 
of the sample to the experimental temperature. It has been stated in 
the Bam® reference that for molybdenum metal and a typical chalcogenii 
glass the ratio of thermal conductivity of metal to that of glass Is 
approximately 500 and this fact makes the infinite heat sink assuap* 
tion an; excellent one. Our glass sample is also sandwiched between 
tt/G. Kot??llic electrodes of noble metal gold. For gold C 

ed.^ CEC Handbook of Chemistry and fhypic's, 61 st edition (CRC 
Boca Batong Florida, 1980- *81) J and soda-lime-silica 
glas ms /** Bef . W* D« .Kingery, Introduction t© Ceramics i Wiley, How 
1967 the ratio of thermal- conductivity of metallic gold to 
®oda-lime-sllica glass Is approximately 190* i’his, therefor®, 
Justifies the boundary condition in section 2.?*4» 
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electric field (E) produces same voltage drop AV - Ew between 
neighbouring layers of neighbouring grains. Here w = is "Uie 

average separation between layers and s^ and d are average interpar- 
ticle separation and diameter of particles respectively. Moreover 
has been assumed that the layers are roughly perpendicular to •ttie 
macroscopic field. Therefore, if t be the ttiickness, V be the 
applied voltage and n be the number of layers. 


E -^ = 
^ - t 


V (V/n) dV 

nTsI+d) s +d s +d 


Therefore applied electric field is uniform throughout the matrix 
and is perpendicular to electrodes. This is one of the assumptions 
(assumption 1, page 25) based on which the formulation has been mac 
Further it should be stated that the computed electric field (B^) 
between the neighbouring metallic particles is assumed to be 
related to the bulk electric field (^) calculated on the basis of 
thickness by the expression (2.39a), page 38. This assumption is 
justified because of the large value of ^^ 1353 / ^etal 
[Ref. 47, page 194], C being the electrical resistivity. 
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It is true thnt wlien a metallic body Is nlnced in a orevlously 

■« ■ , !: 

anlform field between two plates of a large parallel plate capecifcori 
th# field, art und the metallic particle is, distorted. The surface 
remains at constant potential, the field Is perpendicular to equip©- 
tential surface and the force is proportional to surface charge den- | 

i 

gity* But the system which has been investigated is ve^j complicated } 
because of the feet that a large number of metallic particles is i 

distributed in th® glassy phase# Therefore, applying these electro- | 

1 

static principles to an todivldual particle and then extending to the | 
large number of particles present In the system under investigation 
will be an intractable problem to solve. We have therefor# followed I 
an alternative approximate approach which is accepted in literature, j 
£*Eef# B# Abeles, Sheng Ping, M. D. Coutts and Y. Arle, Adv# Phys, | 
2*f, ho? Cl 9755 J7* It; has been stated in literature by Abeles ®t al 
that when an electric field is applied to a granular metal ( metal- I 
dielectric composite ) it is reasonable to assume that the equip©- j 
tential surfaces are approximately perpendicular to the ^plied ; 

macroscopic field and mery metallic particle lies on cm# of th© eipai - 1 
potential surfaces# These surfaces divide the whole granular metal sani 
into a large number of layers so that each particle lying on one of thii 
Issyers la approximately at the same potential. Significant potentl.al 
.'drop, occurs between nelghbourln'f grains only when they are in different 


